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Abstract 

Computation of parallel lines (envelopes) to parabolas, ellipses, and hyperbolas is of importance in 
structure engineering and theory of mechanisms. Homogeneous polynomials that implicitly define paral- 
lel lines for the given offset to a conic are found by computing Grobner bases for an elimination ideal of a 
suitably defined affine variety. Singularity of the lines is discussed and their singular points are explicitly 
found as functions of the offset and the parameters of the conic. Critical values of the offset are linked 
to the maximum curvature of each conic. Application to a finite element analysis is shown. 

Keywords: Affine variety, elimination ideal, Grobner basis, homogeneous polynomial, singularity, 
family of curves, envelope, pitch curve, undercutting, cam surface 
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1 Introduction 

The present paper introduces the Reader to a computation of parallel lines (envelopes) to three nondegenerate 
conies: a parabola, an ellipse, and a hyperbola. In each case, we show how to find a single polynomial 
g G R[x, y] that implicitly generates a suitable elimination ideal of an affine variety containing the parallel 
lines for the given offset r and for the given parameters of the conic. This amounts to a computation of a 
reduced Grobner basis 0], [7], for the ideal for an appropriately chosen total order. We use Maple 
[5] and Singular |7] to accomplish computations. We show an application of the results presented here to a 
finite element analysis. We also make a connection with a design of surface cams in mechanical engineering 

eh, na 

In each of the next three sections we discuss three cases of parallel lines for the specific values of the 
offset < r < r cr i tl r = r cr i t , and r > r cr i t for a parabola (Section^, an ellipse (Section[3J), and a hyperbola 
(Section Then, in each section, we consider a general case of each conic while imposing essentially no 
condition on the offset except that r > 0. In Section [S] we show one possible application of the results to a 
finite clement analysis of a composite material of elliptical shape. In Appendices A and B we collect some of 
the more cumbersome expressions while in Appendix C we show in detail how we found the singular points 
for the parallel lines of the ellipse. 



2 Parallel Lines to a Parabola 

In this section we compute parallel lines to a parabola defined by a polynomial f\ : 

h = ±pyo -4 = (!) 

where \p\ denotes the distance between the focus F — (Q,p) and the vertex V — (0, 0). 1 Polynomial /2 defines 
a circle of radius (offset) r centered at a point (xq, yo) on the parabola f\, 

h = {y - y f + {x ~ x ) 2 - r 2 = (2) 

while a polynomial 

f 3 = 2xp - 2x p + x y - x yo = (3) 

gives a condition that a point P(x, y) on the circle f% lies on a line perpendicular to the parabola f± and 
passing through the point (xo,yo) on the parabola. The family of such circles is called the family of curves 
determined by the polynomial f\ (see 4 ). There are of course two such points for any given point (xq, yo) : 
one on each side of the parabola. All these points P belong to an affine variety V = V(/i,/2,/3) - the 
envelope of the family of circles - and define two parallel lines at the distance r from the parabola. In what 
follows we will compute a reduced Grobner basis for V and by eliminating variables xq, yo we will obtain one 
polynomial g in variables x and y alone (also p, r in a general case) that defines both parallel lines on each 
side of the parabola for the given value of the offset r. Thus, in general, we have /i, /2, f 3 G R[^o, yo, x, y,p, r]. 
However, first we illustrate computation for specific values of the parameters p and r, and later we show 
computation for unspecified parameters. 

1 Without any loss of generality, we assume that the parabola has its vertex at the origin and its focus on the j/-axis. 
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2.1 Parallel lines for specific parameters p and r 

The parallel lines constitute a subvariety of the affine variety V = V(J) where I — (fi, fz, /a). In order to 
find the polynomial g for the second elimination ideal I2 = I H M.[x, y], we first find a Grobner basis G for 
V for the lexicographic order y > x > x > y that contains sixteen polynomials. 2 Then, we reduce G to 
a minimal Grobner basis G m that contains six (or eight, depending on the case) polynomials and, further, 
to a unique reduced Grobner basis G r with three (or five, depending on the case) polynomials only. Finally, 
we eliminate from G r all polynomials that contain variables yo and xo, and keep the single polynomial 
g £ M.[x,y] that gives the parallel lines and generates the second elimination ideal I2 = (g). 

We are also interested when the parallel lines form a smooth (non singular) variety ~V(g) or a C°° 1- 
manifold or a hypersurface defined by g. Since g can be viewed as a C°° function R 2 — > M and V(g) = 
{(x,y) G K 2 \g(x,y) = 0}, it is known that V(g) is a smooth manifold at those points p S V(g) where 
(V<?)(p) 7^ 0. 1101 Theorem 2.6] We will find singular points of V(<?) by finding points p € V(g) where 
(V$)(p)=0. 

Thus, V(g) can be thought of as a union of points where the gradient is zero and those points where the 
gradient is not zero. Connected components of V(g) where the gradient is not zero form the 1-dimcnsional 
manifold of the parallel lines. However, we are also interested in describing a discrete finite set S of singular 
points where the gradient is zero. For a formal definition of a singular point of an affine variety V(<?) see 
page 136]. 

Example 1 (r < r cr it) Let P — \ and r = \. Then, 

V = V(4yo — 3x 2 , 16j/ 2 — 32yyo + 16y 2 + 16x 2 — 32xxo + 16xq — 1, 2x — 2xq + 3xoy — 3xoyo) 
We find that the reduced Grobner basis G r contains only three polynomials including 

g = 331776a; 6 - 42192x 2 + 48400y 2 - 448512y 3 - 84480a; 2 y + 589824y 4 + 28032y 

- 25344a; 4 + 1138176x 2 y 2 - 884736x 2 y 3 - 1105920x 4 y - 5329 + 331776x 4 y 2 (4) 

whose gradient is 

Vg = (-879a; - 1760ya; - 1056a; 3 + 23712xy 2 - 46080a; 3 y - 18432y 3 a; + 13824x 3 y 2 + 20736a; 5 )i 

+ (3025y - 2640a; 2 - 42048y 2 + 876 + 71136a; 2 y - 34560a; 4 - 82944a; 2 y 2 + 73728y 3 + 20736a; 4 y)j . (5) 

Thus, the only real point Si on V(g) where (Vg)(Si) = is Si — (0, j^)- We will refer to this point as a 
virtual singular point for the following reason: V(y) is a disjoint union of {Si} and a non-singular smooth 
subvariety II where the gradient V<? does not vanish. The subvariety II constitutes the set of two parallel 
lines that we have been seeking. We show all in Figure^ 

We have verified by direct computation that as the value of the offset r increases towards some critical 
value r cr a, the virtual singular point approaches the smooth component of the variety V(g), that is, it 
approaches the parallel lines II. 3 As the next example shows, when r = r cr it, then the virtual singular point 
is actually located on the smooth component II of the variety V(g). 

2 Computations have been performed with Groebner package in Maple 8 [5] as well as with a small package written by the 
Authors 0, and the results have been confirmed with Singular 7 using 0. 

3 In the next subsection we will derive a formula for the critical value r cr i t of r for a general parabola. 



3 











i 


















0.8- 






/ g 












0.6- 


















y : 


















0.4- 
\ F 


Sl 
















0.2- 






/ 1 




-1 


-0.8 




-0.4 


-0.2 u 

--___-0.2- 


0.2 


0.4 OS 

X 


O.N 


1 










-0.4- 











Figure 1: Parabola with parallel lines <?, focus F and one virtual singular point Si for P — \ and r = j 

Example 2 (r = r cr u) Let p = -| and r = |. Then, 

V = V(4y - 3a; 2 ,, 9y 2 - 18yy Q + 9y% + 9x 2 - 18a;a;o + 9x 2 - 4, 2x - 2x + 3x y - 3x yo) 
In i/iis case, we /lave found that the reduced Grobner basis G r contains five polynomials including 

g = 768y - 288x 2 - 1728y 3 + 19Ux 2 y 2 - 891x 4 - 2AZQx 4 y 

- 1944x 2 y 3 + 1296y 4 + 729x 4 y 2 + 729x 6 - 256 (6) 

whose gradient is 

Wg = (-32a; + 216a;?; 2 - 198x 3 - 540x 3 y - 216y 3 x + 162xV + 243a; 5 )i 

+ (128 - 864y 2 + 648x 2 y - 405x 4 - 972a;V + 864y 3 + 243x 4 y)j . (7) 

Thus, like in Example^ there is only one real point Si on V(<?) where (Wg)(Si) — : it is Si — (0, |). 
However, unlike in the first example, this time point Si is not disjoint from the rest of the variety V(g) as 
it can be seen in Figured We will refer to this point as a singular point of the variety V(g). The latter can 
still be thought of as union of {Si} and the non-singular smooth subvariety IT where the gradient V g does 
not vanish. We will say that in this case the parallel lines possess one singular point. 

Probably the most interesting case occurs when the value of the offset r > r CT n as shown in the next 
example. In this case, the single singular point of the variety V(<?) splits into three singular points. 
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Figure 2: Parabola with parallel lines g, focus F and one singular point Si for p = | and r = | 

Example 3 (r > r crit ) Lef p = 3 a«rf r = |. T/ien, 

V = V(4y - 3a; 2 ,, Ay 2 - 8yy + Ayl + Ax 2 - %xx + Ax 2 - 9, 2x - 2x + 3x y - 3x Q y a ) 
In this case, we have found that the reduced Grobner basis G r again contains only three polynomials including 

g = 83808y + 52812a; 2 + 16900y 2 - 37248y 3 - 4896x 2 y 2 - 34416a; 4 - 17280x 4 y 

- 13824x 2 y 3 + 92167/ 4 + 5184x 4 j/ 2 + 5184a; 6 - 84681 + 6240x 2 y (8) 

whose gradient is 

\7g = (4401a; - 408xy 2 - 5736a; 3 - 2880a; 3 ?/ - 1152y 3 a; + 864x 3 y 2 + 1296a; 5 + 520yx)i 

+ (10476 + 4225y - 13968?/ 2 - 122Ax 2 y - 2160a; 4 - 5184x 2 y 2 + 4608y 3 + 1296a; 4 y + 780a; 2 )j . (9) 



Thus, unlike in the above two examples, this time there are three singular points: Si = (0, ||) located on the 
y-axis (the symmetry axis of the parabola), and two additional points S2 and S3 



^65 + 36- 123 -27- 12§ 3-123 l 



(10) 



S 3 = - 



^65 + 36-123 - 27-123 3-123 1 



(11) 
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located symmetrically on each side of the y-axis as shown in Figure We will say that in this case the 
parallel lines possess three singular points. A similar figure can be found in 01 page 14 1] 
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Figure 3: Parabola with parallel lines g, focus F and three singular points Si, S2, S3 for p = | and r = | 
2.2 Parallel lines for arbitrary parameters p and r 

In this section we will find a general polynomial g that will generate the second elimination ideal I2 = 
I fl M.[x, y,p, r] where / = (/i,/2,/3) and the generating polynomials are defined in equations |3J), 
and We will also find a formula for the critical value r cr n of r that determines whether V(/2) has one 
or three singular points, and, we will find general formulas for the coordinates of the said singular points for 
arbitrary values of the parameters p and r. 

We begin with a computation of a Grobner basis for the ideal / for the lexicographic order vq > xq > 
x > y > r > p. 

I = (4p£/o - xl,y 2 - 2yy + y\ + x 2 - 2xx Q + x 2 Q - r 2 , 2xp - 2x p + x y - x y ) C R[y ,x ,x,y,r,p] 

Using Maple and Singular PP, [7] we found that a reduced Grobner basis for / consists of fourteen 
homogeneous polynomials of degrees: 7 (two polynomials), 6 (two polynomials), 5 (two polynomials), 4 
(two polynomials), 3 (three polynomials), and 2 (three polynomials). These polynomials are displayed in 
Appendix A. 



G 



Thus, in this general case we observe that I 2 = (g) where g = g\/p is the following homogeneous 
polynomial in R[x, y, r,p] of degree six: 

g = -2pr 2 yx 2 + 8pr 2 y 3 + 8p 2 r 2 y 2 - 32yp 3 r 2 + 16p 4 r 2 - 16yV + 32y 3 p 3 

- I6p 4 y 2 + 3r 2 x 4 + 8p 2 r 4 + 20p 2 r 2 x 2 - y 2 x 4 + lOypx 4 - x 6 - x 4 p 2 + 8py 3 x 2 

- 32x 2 y 2 p 2 + 8x 2 yp 3 - 3r 4 x 2 + 2r 2 x 2 y 2 + r e - r 4 y 2 - 8pr 4 y (12) 

We want to find now the singular points on the variety V = V(g) where the gradient V<? is zero. Thus, V 
can be thought of as union of points where the gradient is zero and those points where the gradient is not 
zero. Connected components of V where the gradient is not zero constitute the 1-dimensional manifold that 
we encountered above. We are interested in the finite set S of the singular points of the variety. We will 
show next that the critical value of r is r crit = \2p\ = 2\p\. Therefore, we compute Vg = Vg(x,y,r,p) : 

V<? = 

< -x(2ypr 2 - 6x 2 r 2 - 20pV + 2x 2 y 2 - 20ypx 2 + 3x 4 + 2x 2 p 2 - 8py 3 + 32y 2 p 2 - 8yp 3 + 3r 4 - 2r 2 y 2 ), 
- x 2 pr 2 + \2y 2 pr 2 + 8yp 2 r 2 - 16p 3 r 2 - 32y 3 p 2 + A8p 3 y 2 - l&yp 4 - yx 4 + hx 4 p + \2x 2 py 2 - Z2x 2 p 2 y 

+ Ax 2 p 3 + 2r 2 yx 2 - r 4 y - Apr 4 , 
r(-2ypx 2 + 8py 3 + 8y 2 p 2 - 32yp 3 + 16p 4 + 3x 4 + 16p 2 r 2 + 20a; 2 p 2 - 6x 2 r 2 + 2x 2 y 2 + 3r 4 

— 2r 2 y 2 — 16ypr 2 ), 

- r 2 yx 2 + Ar 2 y 3 + 8y 2 pr 2 - A8yp 2 r 2 + 32p 3 r 2 - 16y 4 p + 48y 3 p 2 - 32p 3 y 2 + 8pr 4 + 20x 2 pr 2 

+ 5yx 4 - x 4 p + Ay 3 x 2 - 32x 2 py 2 + \2x 2 p 2 y - Ar 4 y > (13) 

In order to find solutions to the vector equation Vg = 0, we need to solve in fact a resulting system of 
four polynomial equations in variables x,y,p,r. The system could be simplified by assuming r > 0, yet to 
efficiently analyze its solutions one can again employ the Grobner basis approach and find a reduced Grobner 
basis for the ideal Is generated by the four equations. Since we want to make sure that the solutions belong to 
V(g), we add equation g = to our system. Let F be a list consisting of the five polynomial equations, that 
is, F = [Vg = 0, g = 0] and let Is = (Vg,g). Thus, we find that the affine variety V(7s) of singular points 
is generated by a reduced Grobner basis consisting of eleven homogeneous polynomials Si,i — 1, ... ,11, in 
y,x,r, and p with total degrees varying from 5 through 17. The third polynomial S3 in the basis has the 
following form: 

s 3 = xr(640/y - 320p 7 + 3r 4 yx 2 - 3r 2 yx 4 - \2pr 2 x 4 + A8x 2 p 3 r 2 - r 6 y + Ur 6 p 

+ 8Ap 3 r 4 - 480/r 2 + yx 6 + I3x 6 p - 2Ax 4 p 3 - 38Ax 2 p b - 96yp 2 r 2 x 2 - 15pr 4 x 2 

- 60yp 2 r 4 + 96yp 4 r 2 - 60yx 4 p 2 - 96x 2 p 4 y) (14) 

It appears that since r and p are non-zero, one possible solution to the equation s 3 = is x = 0. Let's then 
see if the remaining polynomials in F can be solved when x = 0. Upon substituting x = into the list F, 
one obtains 

[0, -(-Apy + r 2 + Ap 2 ){-8py 2 + yr 2 + Apr 2 + Ayp 2 ), r(-Apy + r 2 + Ap 2 )(-2y 2 + 3r 2 - Apy + Ap 2 ), 

A(-y + 2p)(r - y){r + y)(-Apy + r 2 + Ap 2 ), (r - y)(r + y){-Apy + r 2 + Ap 2 ) 2 ] (15) 
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Since r > 0, one possible solution to (|15fl is r = y which gives the following system: 

F 1 = [0, -(r - 2p)\ (r - 2p) 4 r, 0, 0] (16) 

This leads to the critical value of r cr u = 2p. Of course, this makes sense only when p > 0, that is, in general 
we have r cr i t = \2p\ — 2\p\. When p < 0, we can still get the same critical value of r when we select r = — y 
in 115|) which leads to 

F 2 = [0,(r + 2p) 4 r,(r + 2p) 4 r,0,0]. (17) 
Thus, the critical value of r in this case is r cr it — —2p. Combining both cases we get r cr i t = \2p\ = 2\p\. 

Ap 2 -)- T 2 T 2 

Another solution of the system 11511 can be found by setting y = = p H (p ^ 0). This gives 

Ap Ap 

the singular point 

Si = fo, 1 + . AP = (o,p+ when p^0 and r>0 (18) 



Ap J \ Ap / 

found in Examples ^ [21 an( i 121 above. When < r < r crit = 2\p\, then this is the only singular point S\ as 

r 2 

shown in Figure^Jthat is not located on the continuous part of V(g) : It is located at the distance — — above 

A\p\ 

the focus (when p > 0) and below the focus (when p < 0) on the symmetry line of the parabola (the y-axis in 
our case). Then, the parallel line is smooth at the points (0, r) and (0, — r). These points are located on the 
y axis above and below the vertex of the parabola. As r approaches r cr a from the left, then S\ approaches 
one of these points: It approaches the point located on the parallel line on the concave side of the parabola 
(see Figure When r = r cr u then the singular point Si is located on the parallel line on the concave side 
of the parabola (see Figure Q). 

Now we look for those singular points on V(g) for which x ^ 0. These points will be present only when 
t > T cr it. Thus, we now proceed to solve the system of five polynomial equations contained in the list F, 
or, equivalently, the system of eleven polynomial equations Sj = 0, i = 1, . . . , 11, in the Grobner basis for 7g 
under the condition that x ^ 0. Using approach described in more details in Appendix C for the ellipse, in 
the case when r > r cr u one obtains the following two additional singular points 

S*2 = (a, — ) and S3 = ( — a, — ) when p ^ and r > r cr it = 2\p\. (19) 



7/ V 7 

The expressions a, (3 and 7 are as follows 



\ (vr 2 ) i r 2 + 6r 2 » 2 2 a — 3(pr 2 ) 3 pr 2 2s — A(pr 2 ) s p 2 
a =^2 L ^1 i L , (20) 

\Jpr 2 

= pr 2 (22r 6 (pr 2 )i + U52r e p 2 2^ + 7A56p 3 r 4 2% (pr 2 )^ - 6560p 2 r 4 (pr 2 )3 

- 15488p 4 r 4 2^ - 39680p 5 r 2 2^ (pr 2 )^ + 37600p 4 r 2 (pr 2 )i + 7936p 8 23 

- 15872p 7 2i(pr 2 )^ +45760p 6 r 2 2^ - 33r 8 2^ - 5376p 6 (pr 2 )i), (21) 

7 = 2(pr 2 )^(mA0p 6 r 2 {pr 2 )i - 3968p 7 23 ( pr 2 )i + 1984p 7 r 2 2§ 

- 3920rV(pr 2 )3 - 9920p 5 r 2 23 ( pr 2 )i + 4960p 5 r 4 2§ + A8Ar 6 p 2 (pr 2 )i 

+ imAp 3 r 4 2z (pr 2 )i - 932p 3 r 6 2s - llr 8 (pr 2 )s). (22) 

Thus, in summary, we have these three cases: 
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(i) One (virtual) singular point Si = (0,p + ^) when r < r cr u — 2\p\. In particular, when P = \ and 
r = |, we get from <|18ll that Si = (0, j^) : That is the only (discrete) singular point of V(g) found in 
Example n Formula (|12|) then gives 



5329 3025 , 73 55 , 73 o 293 , 

a = y y A yx A y H x 

y 331776 20736 y 864 y 216 y 54 y 2304 

11 4 247 2 2 16 4 10 4 8 3 2 6 2 4 

144 72 y 9 y 3 y 3 y y 

which is, up to the form and sign, the same as polynomial 

(ii) One (true) singular point Si = (0,p+ ^) when r = r cr ;t = 2|p|. In particular, when P = \ and r = |, 
we get from ()18ll that Si = (0, |) : That is the only singular point of located on the continuous part of 
V(y) found in Example [21 Formula (|12l) then gives 

256 256 64 o 32 , 11 4 8 , 2 16 4 10 4 8 „ 2 fi 2 4 

q = y A u H a; A x x y y A yx + -y J x — x — y x 

y 729 243 y 27 y 81 9 3 y 9 y 3 y 3 y y 

which is, up to the form and sign, the same as polynomial l{6*|). 

(hi) Three singular points: Si = (0,p+ S 2 = (a, and S3 = (-a, |) when r > r c „ t = 2|p|. In 
particular, when P — \ and r = |, we get from (I18H that Si = (0, ||) while from l|19|l we get the two 
additional singular points S 2 and S3 displayed in (|1U|) and (|llf> found in Example 3. Formula l|12f) then 
gives 

9409 4225 2 97 65 2 194 . 163 2 239 4 

q = y y yx A y x A x 

576 1296 y 6 y 54 y 27 y 16 36 

IT 9 9 16 a 10 4 809 ^ 94 

H ari/ y 4 H yx 4 + -y 6 x z - x - y 2 x 4 

18 y 9 y 3 y 3 y y 

which is, up to the form and sign, the same as polynomial (|8"|). 

As a final comment, let's notice that for the parabola 4p(y — k) = (x — h) 2 the critical value of the offset 

Tcrit 2jp| pmin 
^raax 

is the reciprocal of n max , the maximum curvature of the parabola at its vertex, or, equivalently, it equals 
Pmin, the minimum radius of the osculating circle. Furthermore, the quantity 2\p\ is known as the semi-latus 
rectum and represents the distance from the parabola focus to the parabola measured along a line parallel 
to the parabola directrix. ^3] 

3 Parallel Lines to an Ellipse 

In this section we compute various parallel lines to an ellipse defined by a polynomial 

fi = b 2 ylAra 2 xl~ a 2 b 2 =0 (23) 

Recall that when /3(t) is a regular curve in R 3 then k= \$ X /3|/|/3| 3 . EH page 46]. 



9 



where we assume that a > b and the center of the ellipse is at the origin. Polynomial fa defines a circle of 
radius (offset) r centered at a point P = (xo, yo) on the ellipse, 



fa = (y - yo) 2 + (x- x ) 2 - r 2 = 



(24) 



while polynomial fa 



fa = b 2 (x - x )y - a 2 (y - y )x a = 



(25) 



gives a condition that a point P(x, y) on the circle fa lies on a line perpendicular to the ellipse fa and passing 
through the point (xo, yo) on the ellipse. All these points P belong to the affine variety V = V(/i, fa, fa) 
- the envelope - and, like for the parabola, define two parallel lines at the distance r from the ellipse. We 
proceed like in Section [21 We will compute a reduced Grobner basis for V and by eliminating variables 
Xo, yo we will obtain one polynomial g in variables x and y alone (also a,b,r in a general case) that defines 
both parallel lines on each side of the ellipse for the given value of the offset r. Thus, in general, we have 
/li fa, fa € R[?/0) ^Oj x, y, a, b, r\. However, first we illustrate computations for specific values of the parameters 
a, b and r, and later we show computation for unspecified parameters. 

3.1 Parallel lines for specific parameters a, b, r 

The parallel lines constitute a subvariety of the affine variety V = V(I) where, like before, I(fi, fa, fa). 
In order to find the polynomial g for the second elimination ideal fa = I H R[x, y], we first find a Grobner 
basis G for V for the lexicographic order yo > xq > x > y and reduce it eventually to a unique reduced 
basis G r of nine (or eleven, depending on the case) polynomials only. Finally, we eliminate from that basis 
all polynomials that contain variables yo and xo, and keep the single polynomial g £ R[x,y] that gives the 
parallel lines and generates the second elimination ideal fa = (g). Wc will also proceed to find all singular 
points p of V(<7) where (Vg)(p) = 0. Thus, let's consider three cases when < r < r crit , r = r crit , and 
f > TcHt- In the next section we will consider a general case for arbitrary values of a, b (a > b > 0) and 
r (r > 0). 

Example 4 (r < r C rit) Let a = 3, b = |, and r = h. Then, 



g = 256x 8 + 2080x 6 - 41685x 2 + 44100 - 25356y 2 + 5353y 4 - 4751x 4 - 488y 6 - 3360y 4 x 2 

- 4680y 2 x 4 + 528y 4 x 4 + 640x 6 y 2 + 160y 6 x 2 + 16y 8 + 21410x 2 y 2 (26) 



V = \(yl + 4x^ - 9, 4y 2 - 8yy + 4y% + 4x 2 - 8xx + 4x^ - 1, y x + 3y x ~ 4x y) 



We find that the reduced Grobner basis G r contains nine polynomials including 



whose gradient is 



Vg = x(1024x 6 + 6240x 4 - 41685 - 9502x 2 - 3360y 4 - 9360x 2 y 2 + 1056y 4 x 2 

+ 1920y 2 x 4 + 160y 6 + 21410y 2 )i 
+ y(-12678 + 5353y 2 - 732y 4 - 3360x 2 y 2 - 2340x 4 + 528y 2 x 4 + 320x 6 



(27) 



+ 240y 4 x 2 + 32y 6 + 10705x 2 )j . 



There are exactly two solutions to the equation (X7g)(p) — on *V(g) : 



Si = (0, V6) and S 2 = (0, - V6) 
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We will refer to these points as virtual singular points for the following reason: V(g) is a disjoint union 
of {Si, S2} and a non-singular smooth subvariety H where the gradient Vg does not vanish. The subvariety 
II constitutes the set of two parallel lines that we have been seeking. Note that the foci F\ — (0, ^y^) and 
F2 = (0, — ^rp) are different from the singular points. We show all in Figure ^| 



4i 




-4 J 



Figure 4: Ellipse with parallel lines g, foci F%, F2 and two singular points Si, S2 for a = 3, b = | and r = i 

We have verified by direct computation that as the value of the offset r increases towards some critical 
value r cr it, the virtual singular points approach the smooth component of the variety V(^), that is, they 
approach the parallel lines II. 5 As the next example shows, when r = r cr it, then the virtual singular points 
are actually located on the smooth component II of the variety V(g). 

Example 5 (r = r cr n) Let a = 3, b = |, and r = §. Then, 

V = V(y 2 + 4a; 2 _ 9, i6 y 2 - 32 TOo + 16y 2 + 16x 2 - 32a;a; + 16x 2 - 9, y Q x + 3y Q x - Ax y) 

We find that the reduced Grobner basis G r contains eleven polynomials including 

g = -29673216a; 4 + 19035648y 4 + 95178240a; 2 y 2 - 1916928/ - 14376960yV 

- 21012480?/ 2 a: 4 + 7372800a: 6 + 65536y 8 + 2162688yV + 655360j/V 

+ 2621440a;V + 1048576a; 8 - 79361856?/ 2 + 119574225 - 198404640a; 2 (28) 

5 In the next subsection we will derive a formula for the critical value r c rit °f r f° r a general ellipse. 
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(29) 



whose gradient is 

Vg = a;(-1854576x 2 + 2974320?/ 2 - 449280y 4 - 1313280x 2 y 2 + 691200a; 4 + 135168yV 

+ 20480y 6 + 245760yV + 131072a; 6 - 6200145)i 
+ y(594864y 2 + 1487160a; 2 - 89856y 4 - 44928CteV - 328320a; 4 + 4096y 6 

+ 67584yV + 30720y 4 a; 2 + 40960a; 6 - 1240029)j . 

There are exactly two solutions to the equation (Vy)(p) = on V(<?) : 

( 9 \ ( 9 

Si = I 0, -J and S , 2=(0,— - 

These points are the singular points of the variety V(g) : At every other point of this variety the gradient 
Vg does not vanish. We will collect these points as before into a non-singular subvariety II. The variety 
II constitutes the set of two parallel lines that we have been seeking. Note that the foci F\ = (0, ^y^) and 



F2 = (0, 



3yg 



) are still different from the singular points. We show all in Figure 




-4- 



Figure 5: Ellipse with parallel lines g, foci F±, F2 and two singular points Si, S2 for a = 3, b = | and r — 

3 

— ' crit — 

are located on the continuous component of II. 



As we can see from the last figure, when r = r cr u = | for the given values of a and b, the singular points 
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Probably the most interesting parallel lines we obtain when r > r cr i t as we show in the next example. 
There, we will find six singular points as each of the two singular points found in Example [5] splits into three 
points. 



V = V(y 2 + 4a; 2 - 9, 9y 2 - 18yy + 9y% + 9x 2 - 18xx + 9a; 2 - 16, y x + 3y x - Ax y) 
We find that again the reduced Grobner basis G r contains nine polynomials including 

g = 23461560a; 2 ?/ 2 + 1221025 + 186624a; 8 - 2984040y 4 a; 2 - 5015520yV 

- 284472y 6 + 384912?/V - 43658160a; 2 + 518400a; 6 + 466560a; 6 ?/ 2 

+ 116640y 6 a; 2 + 11664y 8 - 2228394?/ 2 - 10769184a; 4 + 1344177?/ 4 (30) 

whose gradient is 

V.</ = a;(977565?/ 2 + 31104a; 6 - 124335?/ 4 - 417960a; 2 ?/ 2 + 32076?/ 4 a; 2 - 1819090 + 64800a; 4 



There are exactly six solutions to the equation (S7g)(p) = on V((/), or, in another words, six singular points 
of the variety: 



Example 6 (r > r cr jt) Let a — 3, b — |, and r = |. Then, 



+ 58320?/ 2 a; 4 + 4860?/ 6 - 897432a; 2 )i 
+ ?/(3910260a; 2 - 994680x 2 y 2 - 835920a; 4 - 142236?/ 4 + 128304?/ 2 a; 4 + 77760a; 6 

+ 58320?/ 4 a; 2 + 7776y 6 - 371399 + 448059?/ 2 )j . 



(31) 





We display these six points in Figure^ 
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Figure 6: Ellipse with parallel lines <?, foci Fi, F2 and six singular points Si, i = 1, . . . , 6, for a = 3, b = | 
and r = | 

3.2 Parallel lines for arbitrary parameters a, b, r 

In this section we will find a general polynomial g that will generate the second elimination ideal I2 = 
I D M.[x, y, a, b, r] where I = (fi, f2, /a) and the generating polynomials are defined in equations (|23|l . 12411 . 
and l|25|l . We will also find a formula for the critical value r cr i t of r that determines whether V(J2) has two 
or six singular points, and, we will find general formulas for the coordinates of the said singular points for 
arbitrary values of the parameters a, b and r. 

We begin with a computation of a Grobner basis for the ideal / for the lexicographic order yo > xq > 
x>y>a>b>r. 

I = (b 2 yl + a 2 xl - a 2 b 2 ,y 2 - 2yy + y 2 + x 2 - 2xx + x 2 - r 2 ,b 2 y Q x - b 2 y x Q - a 2 x y + a 2 x y ) 

We found that a reduced Grobner basis for / consists of fifteen homogeneous polynomials of degrees: 16 
(two polynomials), 14 (one polynomial), 12 (one polynomial), 11 (two polynomials), 10 (two polynomials), 
9 (one polynomial), 8 (one polynomial), 7 (two polynomials), 5 (one polynomial), 4 (one polynomial), and 2 
(one polynomial). These polynomials are displayed in Appendix B. 

Thus, in this general case we observe that I2 = (g) where g = gi/(a 2 b 2 ) is the following homogeneous 
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polynomial in R[x, y, a, b, r] of degree twelve 6 : 



g = 6b 4 x 4 a 4 - 4y 6 b 4 a 2 - 6 a 6 x 4 b 2 - 46 6 x 2 a 4 + 6 b 4 x 2 a 6 + a 4 b 8 - 2 y 2 b 8 a 2 

+ 6y 2 b 6 x 2 a 2 + 6 a 6 y 2 x 2 b 2 - 10 b 4 x 2 a 4 y 2 + y 8 b 4 + b 8 r 4 + b 4 r 8 - 2 b 6 r 6 

- 2 b 4 y 2 x 4 r 2 - 6 y 4 b 4 r 2 + 6 x 2 6 4 r 4 y 2 ~ 4 6 4 y 6 r 2 + 2 6 4 y 6 x 2 + y 4 b 4 x 4 

+ b 4 x 4 r 4 -2b 4 r 6 x 2 -2b 6 x 2 r 4 + 6y 4 b 4 r 4 -2b 8 y 2 r 2 + 6b 6 r 4 y 2 -4b 4 r 6 y 2 
+ 6 a 4 y 2 6 6 + 6 a 6 y 2 6 2 r 2 + 2 a 4 6 2 x 4 y 2 - 4 a 4 b 2 + 4 a 6 x 2 b 2 r 2 - 8 6 4 x 2 a 4 r 2 
+ 10 a 4 x 4 b 2 r 2 - 6 b 4 r 2 a 2 x 2 y 2 - 6 a 4 x 2 b 2 r 2 y 2 - 8 a 4 x 2 6 2 r 4 - 6 6 4 r 2 a 2 x 4 
+ Ab 4 r 4 a 2 x 2 ~& y 4 a 4 b 2 x 2 - 6 6 4 x 4 a 2 y 2 + 6 6 6 x 2 r 2 a 2 + 2 6 4 x 2 y 4 a 2 
+ 10b 4 y 4 a 2 r 2 +4b 6 y 2 r 2 a 2 - 8 6 4 y 2 r 4 a 2 - 8 6 4 y 2 a 4 r 2 + 2 r 6 a 2 6 4 + a 8 6 4 

- 6 a 4 y 4 r 2 b 2 + 4 a 4 y 2 r 4 b 2 - 2 a 8 6 2 x 2 + x 4 a 4 y 4 + 4 z 2 6 6 y 2 r 2 - 6 y 4 b 6 a 2 (32) 

- 4 a 6 y 2 6 4 - 2 y 4 a 4 x 2 r 2 + 2 6 2 r 2 a 2 x 4 y 2 - 2 r 8 a 2 6 2 + r 8 a 4 + a 8 + a 4 x 8 
-4a 4 x 6 r 2 + 6a 4 x 4 r 4 + 2 b 2 r 2 a 2 x 6 - 6 6 2 r 4 a 2 x 4 - 2 a 6 r 6 + 2 a 6 x 6 + a 8 r 4 

,o6i,2 22 o6 4 2.c6i.2 2 2 o822, 444 c 4j2 42 

+ 2y b a x — 2 r a y + 6 r o y a — 2 a .t r + y a r — b y b r a 
-2a 6 y 2 r 4 + 4x 2 a 6 y 2 r 2 + 6 a: 2 a 4 y 2 r 4 + 2 y 4 b 2 x 2 a 2 r 2 - 10 x 2 b 2 y 2 r 4 a 2 
+ 2y 6 b 2 a 2 r 2 + 6 r 6 b 2 x 2 a 2 + 4 b 2 x 4 y 4 a 2 + 2 a 4 x 6 y 2 - 6 a 6 x 4 r 2 - 4 r 6 a 4 x 2 
+ 26 2 a ; 6 a 2 y 2 + 6r 4 a 6 2 ; 2 -6a 4 a; 4 r 2 y 2 + 26 6 y 6 -2 a 8 6 2 r 2 -2 a ; 2 y 4 6 6 
+ 6a 4 y 4 6 4 + y 4 6 8 -6y 4 6 6 r 2 -2 a; 4 a 6 y 2 -2a 6 6 6 + 26 6 r 4 a 2 -26 8 r 2 a 2 
+ 2 b 4 a 6 r 2 - 6 b 4 a 4 r 4 + 2 r 2 a 4 b 6 + 2 r 4 a e b 2 + 2 r e a 4 b 2 

We want to find now those points on the variety V = V(g) where the gradient Vg is zero. Thus, V can 
be thought of as union of points where the gradient is zero and those points where the gradient is not zero. 
Connected components of V where the gradient is not zero constitute the 1-dimensional manifold. We are 
interested in a discrete finite set S of points where the gradient is zero. These points will give the singular 
points of the variety V. We will show next that the critical value of r is r cr i t = Therefore, we compute 
the gradient Vg = Vg(x, y, a, b, r) =< h\, h 2 , h 3 , hi, h 5 > where 

/ il =x(2 a 4 6 6 -2 a 4 6 2 a ; 2 y 2 -66 4 a; 2 a 4 + 56 4 y 2 a 4 + 6 6 y 4 + 66 4 y 2 a 2 a ; 2 + 6 6 r 4 

- 2 b 6 y 2 r 2 - a 8 x 2 - 2 a 6 y 2 r 2 + a 4 y 4 r 2 - 3 a 4 y 2 r 4 + 2 a 4 r e + a 8 r 2 - 3 a 6 r 4 

- b 4 r 4 x 2 +b 4 r 6 -A x 2 b 2 a 2 y 4 - 3 x 4 b 2 a 2 y 2 - x 2 b 4 y 4 - a 2 y 6 b 2 - 3 x 4 a 4 y 2 

+ 6a 4 x 2 r 2 y 2 - b 2 y 4 a 2 r 2 + 5 b 2 y 2 r 4 a 2 - 3 b 2 a 2 x 4 r 2 + 6 b 2 a 2 x 2 r 4 + 6 a 4 x 4 r 2 

+ 6a 6 x 2 r 2 -6a 4 x 2 r 4 + 2 b 4 y 2 x 2 r 2 - 3 b 4 r 4 y 2 + 3 y 4 b 4 r 2 - 3 b 2 r 6 a 2 (33) 

-2 b 2 a 2 x 2 r 2 y 2 + 3 y 4 a 4 b 2 + a 8 b 2 - 2 a 6 b 2 r 2 + 4 a 4 b 2 r 4 + 4 a 4 b 4 r 2 - 3 a 6 x 4 

- b 4 y e + 3 b 4 y 2 a 2 r 2 - 10 a 4 b 2 x 2 r 2 + 3 a 4 b 2 y 2 r 2 + 6 b 4 x 2 a 2 r 2 + 6 x 2 a e b 2 
-3a 6 b 2 y 2 -3b e y 2 a 2 - y 4 b 4 a 2 - 3 b e r 2 a 2 - 2 b 4 r 4 a 2 + 2 x 2 y 2 a 6 + 6 x 4 a 4 b 2 

-3a 6 b 4 -a 4 y 4 a; 2 -2a ; 6 a 4 ) 

6 We assume a > b > 0. 
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h 2 =y(3a 2 b 4 x 4 -3aH 6 + 6a 4 b 2 x 2 y 2 + 5b 4 x 2 a 4 - 6 b 4 y 2 a 4 + 2 b 6 y 2 x 2 - 3 b 6 y 4 

- 2 6 4 y 2 a 2 x 2 - 3 b 6 r 4 + 6 6 6 y 2 r 2 - a 4 y 2 r 4 + a 4 r 6 + a 6 r 4 + a; 4 fe 4 r 2 - 3 6 4 r 4 x 2 
-2b 6 x 2 r 2 + b s r 2 + 2b 4 r 6 -3 x 2 b 2 a 2 y 4 - 4 a: 4 6 2 a 2 y 2 - x 6 b 2 a 2 - 3 x 2 6 4 y 4 

- x 4 b 4 y 2 ~x 4 a 4 y 2 + 2 a 4 x 2 r 2 y 2 - 3 b 2 y 4 a 2 r 2 + 6 6 2 y 2 r 4 a 2 - 6 2 a 2 x 4 r 2 

+ 5b 2 a 2 x 2 r 4 + 3 a 4 x 4 r 2 - 2 a 6 x 2 r 2 - 3 a 4 a; 2 r 4 + 6 6 4 y 2 x 2 r 2 ~ 6 6 4 r 4 y 2 (34) 

+ 6 y 4 b 4 r 2 - 3 b 2 r 6 a 2 - 2 b 2 a 2 x 2 r 2 y 2 + b s a 2 - 3 a 6 b 2 r 2 - 2 a 4 6 2 r 4 

+ 4a 4 6 4 r 2 + a 6 a; 4 -6 8 y 2 -26 4 y 6 -106 4 2 ; 2 a 2 r 2 + 3a 4 6 2 a; 2 r 2 + 6a 4 6 2 ?/ 2 r 2 

+ 3 b 4 x 2 a 2 r 2 - 3 x 2 a e b 2 - 3 x 2 b 6 a 2 + 6 6 6 y 2 a 2 + 6 y 4 6 4 a 2 - 2 6 6 r 2 a 2 

+ 4 6 4 r 4 a 2 - x 4 a 4 b 2 + 2 a 6 6 4 - x 6 a 4 ) 

/i 3 = 6a 2 6 4 2 ; 4 + 3 a 4 6 6 -9 a 4 6 2 a ; 2 y 2 - 96 4 a; 2 a 4 + 66 4 y 2 a 4 -36 6 y 2 a ; 2 + 36 6 y 4 

+ 106 4 y 2 a 2 x 2 -6 6 r 4 -26Vr 2 + 3a 4 y 2 r 4 + 3a 4 r 6 -2a 6 r 4 + 3x 4 6 4 r 2 
-2b 4 r 4 x 2 -3b 6 x 2 r 2 + b s r 2 - b 4 r 6 + 6 x 2 b 2 a 2 y 4 ~ 2 x 4 b 2 a 2 y 2 + 4 x 6 6 2 a 2 

- x 2 b 4 y 4 + 3 x 4 b 4 y 2 -a 2 r 8 -3 r 6 x 2 b 2 - 3 r 6 b 2 y 2 + 5 x 2 6 2 y 2 r 4 - x 4 6 2 r 2 y 2 

2 j.2 4 2 28 6 j.2 2 <-, 4 j2 4 e 2 2 2 4 6 1,2 2 , o 4 i2 4 

- x o y r —ax— x b y — 2 x of - 6 a x y r — y o r + 3 y b r 
+ 6 a 2 x 4 r 2 y 2 + 2 a 2 y 4 x 2 r 2 - 6 2 x 2 y 6 - x 6 fe 2 r 2 + 3 x 4 b 2 r 4 + 3 x 4 a 4 y 2 

- 6 a 4 x 2 r 2 y 2 + 6 fe 2 y 4 a 2 r 2 - 4 fe 2 y 2 r 4 a 2 - 10 6 2 a 2 x 4 r 2 + 8 fe 2 a 2 x 2 r 4 (35) 
+ 9a 4 x 4 r 2 +4a 6 x 2 r 2 - 9a 4 x 2 r 4 + 3 b 4 y 2 x 2 r 2 + Ab 4 r 4 y 2 -hy 4 b 4 r 2 

-2b 2 r 6 a 2 + 6b 2 a 2 x 2 r 2 y 2 -b 8 a 2 + 4a 6 b 2 r 2 - 3 a 4 b 2 r 4 - 3 a 4 b 4 r 2 - 2 a 6 x 4 

+ b s y 2 + 2 b 4 y 6 + 8 b 4 y 2 a 2 r 2 - 6 a 4 b 2 x 2 r 2 - 9 a 4 b 2 y 2 r 2 + 8 b 4 x 2 a 2 r 2 

+ 4 x 2 a 6 b 2 + 4 x 2 b 6 a 2 - 6 6 6 y 2 a 2 - 6 y 4 fe 4 a 2 - 2 b 6 r 2 a 2 + 6 6 4 r 4 a 2 

+ 9 x 4 a 4 b 2 - 2 a 6 b 4 - 3 x 6 a 4 + 4 a 2 x 6 r 2 + 4 a 2 r 6 x 2 - a 2 y 4 x 4 - 6 a 2 x 4 r 4 

- 2 a 2 x 6 y 2 - a 2 y 4 r 4 + 2 a 2 r e y 2 + r 8 b 2 
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h 4 = 2a i b 6 -WaH 2 x 2 y 2 -6b i x 2 a 4 + 9b 4 y 2 a 4 + 2 b 6 y 4 + 9 b 4 y 2 a 2 x 2 + 2b 6 r 4 

-4b 6 y 2 r 2 -a 8 x 2 + 3 a 6 y 2 r 2 - 3 a 4 y 4 r 2 + 2 a 4 y 2 r 4 + a 4 r 6 - a 8 r 2 + a 6 r 4 
-3b 4 r 4 x 2 -3b 4 r 6 + 2x 2 b 2 a 2 y 4 - 6 x 4 b 2 a 2 y 2 - 3 x 2 6 4 y 4 - 4 a 2 y 6 6 2 - a 2 r 8 

n SJi2 /i6i.2 2ic2i,2 24 o4i2 22 c ™2 j.2 „ 4 2 , 4 i.2 4 

— 2 r x o — 4 r b y + 6 x b y r — 2 or o r y — 6 x b y r + x b y 

r-2224 /l 6 1.2 „2 , c 4 j.2 „4 i 2422, 2 4 22,r,7.2 2 6 

— 5 a x y r — 4 y o r + 6 y o r + a x r y + a y x r + 2 o a; y 

i 4l2 4, 442 o4222|-ir>J.2 422 o 1.2 242 e 1.2 242 

+ at o r + x a y — 3 a x r y + 10 o y a r — 8 b y r a — 6 o a x r 

+ 4 6 2 a 2 x 2 r 4 + 5 a 4 x 4 r 2 + 2 a 6 x 2 r 2 - 4 a 4 x 2 r 4 + 6 6 4 y 2 x 2 r 2 + 9 6 4 r 4 y 2 (36) 

— 9 y 4 6 4 r 2 + 2 b 2 r 6 a 2 - 6 b 2 a 2 x 2 r 2 y 2 + 6 y 4 a 4 6 2 + a 8 6 2 + 2 a 6 fe 2 r 2 
-6a 4 6 2 r 4 + 3a 4 6 4 r 2 -3a 6 x 4 + 36 4 y 6 + 66 4 y 2 a 2 r 2 -8a 4 6 2 x 2 r 2 
-8a 4 6 2 y 2 r 2 +9b 4 x 2 a 2 r 2 + b 2 y 8 + a 2 y e r 2 + x 2 a 2 y 6 + 6x 2 a 6 6 2 -4a 6 b 2 y 2 

— 4 6 6 y 2 a 2 - 9 y 4 b 4 a 2 - 4 6 6 r 2 a 2 + 3 6 4 r 4 a 2 + 3 x 2 y 2 a 6 + 6 x 4 a 4 b 2 - 3 a 6 6 4 

o442 r,64, 262io262,r,244 o244, 262 

— 3 a y x — 2 x a + a x r + 3 a r x + 2 a y x — 3 a x r + a x y 

o244io262, 8 t2 

— 3 a y r + 3 a r y + r b 



(37) 



h 5 =3a 2 b 4 x 4 + a 4 b 6 ~3 a 4 b 2 x 2 y 2 - 4 b 4 x 2 a 4 - 4 6 4 y 2 a 4 + 2 6 6 y 2 x 2 - 3 b 6 y 4 

-3b 4 y 2 a 2 x 2 -3b e r 4 + 6 b 6 y 2 r 2 - a 8 x 2 - 2 a 6 y 2 r 2 + a 4 y 4 r 2 - 3 a 4 y 2 r 4 
+ 2 a 4 r 6 + a 8 r 2 - 3 a 6 r 4 + x 4 b 4 r 2 - 3 b 4 r 4 x 2 - 2 b 6 x 2 r 2 + b 8 r 2 + 2 b 4 r 6 

i 2j,2 24, 4 1,2 22, 6 l2 2 o2j,4 4 4 l4 2i 26i,2 o442 

+ x o a y + x b a y + x b a — 3 x o y — x b y + a y b — 3 x a y 

, e 4222 ^ 7.2 422,^1.2 242 c i2 242,riJ,2 224 

+ 6a x r y —66 y a r +9b y r a -do a x r + 9 o a x r 

ie442,c622 a 4 _2 4 , a J.4 2 „2 2 oj.4 42, c 4i,4 2 

+ 6 a x r + 6 a x r — 6 a x r + 6 b y x r — 6 o r y + 6 y o r 

,1 1.2 62 -inj.2 2222 o44l2 8 j.2 l8 2 , r> 6 1.2 „2 i o 4 j.2 „4 

— 4 o r a — 10 o a x r y — 3 y a o — a o — o a + 2 a b r + 3 a b r 

- 6a 4 6 4 r 2 - 3a 6 x 4 - b 8 y 2 - 26 4 y 6 - 86 4 y 2 a 2 r 2 - 8a 4 6 2 x 2 r 2 + 4a 4 b 2 y 2 r 2 
+ 4 & 4 x 2 a 2 r 2 + 2 x 2 a 6 b 2 + 3 x 2 b 6 a 2 + 3 a 6 b 2 y 2 + 2 b 6 y 2 a 2 + 5 y 4 6 4 a 2 

+ 2b 6 r 2 a 2 + 3b 4 r 4 a 2 + 2 x 2 y 2 a 6 + 5 x 4 a 4 6 2 + a 6 b 4 - a 4 y 4 x 2 - 2 x 6 a 4 

are polynomials in K[a, y, a, 6, r]. Continuing our approach from Section [21 in order to find solutions to the 
vector equation Vy = in V(y), we need to solve a system of six polynomial equations hi = 0, i = 1, . . . , 5, 
and g = in the variables x,y,a,b,r under the assumption that a > b > r > 0. Let Is — (Vy, y) = 
(/ii, /12, /14, <?) and let F be a list containing the six polynomial equations: 

F = [hi = 0, h 2 = 0,h 3 = 0, hi =0,h 5 = 0, g = 0] (38) 

In order to solve system (I38|l . we can either employ the Grobner basis approach again (see Appendix C) or 
simply use solve command in Maple. Thus we find that the affme variety V(/s) consists of the following 
eight singular points of which six are real and two are complex: 



a _ LVE^EnEEn , a _ ,_v2EZSFZS , (39) 
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«■=-,- , , .... . -a .- ,r\, : ' ■ («) 



(-a 2 + 6 2 ) (a6r)f ' 


(a6r)s (-a 2 + & 2 ) y 




V(-a 2 + 6 2 )e 


(-a 2 + 6 2 ) (a6r)§ ' 


(abr)i (-a 2 + 6 2 ) 






(-a 2 + b 2 ) (abr)i ' 


(a6r)i (-a 2 + 6 2 ) y 


/ 


V(-a 2 + 6 2 )e > 



S V = ( » — — • ^ / ', • ) • (43) 

I • ( 44 ) 
(—a 2 + 6 2 ) (a 6 r) 3 (a 6 r) 3 (— a 2 + & 2 ) y 



where 



<5 = -3 b 2 a 2 r 2 + r 2 a 2 (abr)t +3{abr)^ b 3 ar ~ b A (abr)T> , (45) 
e= -a 4 (a6r)i + 3 (a br)H r a 3 - 3 b 2 a 2 r 2 + (a br)H 2 r 2 . (46) 
Notice that 5 becomes zero in three cases: 

1. When r = - this is not possible since we expect r > 0, 

,2 

2. When r = — — - this is not possible as r > 0, 

3. When r = ^. 

a 

,2 

Only this last value is possible thus we set r cr n = — . When r = r cr i t , then under the assumptions that 
a > b > r > 0, we essentially obtain two real singular points: 

/ a 2 - b 2 \ ( a 2 - b 2 \ 
Si = S 5 = S 6 = f 0, — - — J and S 2 = 5V = S 8 = f 0, — J , (47) 

while the complex points S3 and S4 become: 

_ / V^TF(a 2 - b 2 ) . \ / V^T¥(a 2 -b 2 ). \ , AQ ^ 
S 3 = '-%, I and Si = I ^, I . (48) 

For the ellipse with a > 6, the ratio / = ^- is again called a semi-latus rectum. It is the distance between 
a focus and the ellipse itself measured along a line perpendicular to the major axis. The value of I is also 
the reciprocal of the maximum curvature n max of the ellipse, or, it is equal to the radius of the smallest 
osculating circle. [Tl] 

To summarize, the variety V(is) contains the following real singular points: 

1. Two virtual points S\ and S2 shown in (|39|) when < r < r cr u. In particular, when a = 3, b = |, and 
r = |, upon substituting these values into ll-'SL'l we obtain the parallel lines (|26() derived in Example 0] 
Furthermore, when these values are substituted into l|39l) . then we obtain the two virtual singular 
points found in the same example. 
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2. Two points S\ and S*2 shown in (|47|l when r — r cr i t . In particular, we can recover results of Example^] 
when we substitute values a = 3, b = I, and r = | into the general equation ilMl'l) and the singular 
points lET7|) . 

3. Six points Si, S2, S5, Sq, Sr, S$ shown in JSHl, JUJ, 621, (E3), and JUJ when r > r crlt . In particular we 
can recover in the same manner six points and the equation for the parallel lines found in Example El 
when we set a = 3, b = |, and r = §. These six singular points are displayed in Figure El 

4 Parallel Lines to a Hyperbola 

Since a parabola, an ellipse, and a hyperbola belong to the same class of nondegenerate conic 

C : (X 2 + Y 2 - Z 2 = 0) 

in P 2 (M), (cf. we expect essentially similar results for a hyperbola. Thus, our approach will be similar to 
the one for the ellipse from Section|31 In particular, there is only one sign change needed in the polynomials 
fa and fa shown in i12.'j[| and (|25|l whereas the polynomial fa remains the same. Thus, we have now 

fa = b 2 y 2 - a 2 x 2 - a 2 b 2 = (49) 

where we assume that a > b and the center of the hyperbola is at the origin. Polynomial fa defines a circle 
of radius (offset) r centered at a point P = (xq, yo) on the hyperbola, 

fa = (V - yof + (x- x ) 2 -r 2 =Q (50) 

while polynomial fa 

fa = b 2 (x - x a )y + a 2 (y - y )x = (51) 

gives a condition that a point P(x, y) on the circle fa lies on a line perpendicular to the hyperbola fa 
and passing through the point (xo,yo) on the hyperbola. All these points P belong to the affine variety 
V = V(/i,/2,/3) - the hyperbola envelope - and, like for the parabola in Section |21 and the ellipse in 
Section |3 define two parallel lines at the distance r from the hyperbola. We proceed like in Section |3J We 
will compute a reduced Grobner basis for V and by eliminating variables xo , yo we will obtain one polynomial 
g in variables x and y alone (also a, b, r in the general case) that defines both parallel lines on each side of 
the hyperbola for the given value of the offset r. Thus, in general, we have fa, fa, fa G M.[yo, xo, x,y,a,b,r]. 
However, first we illustrate computations for specific values of the parameters a, b and r, and later we show 
computation for unspecified parameters. 

4.1 Parallel lines for specific parameters a,b,r 

The parallel lines constitute a subvariety of the affine variety V = V(I) where, like before, I(fi, fa, /s). In 
order to find the polynomial g for the second elimination ideal fa = lC\M.[x, y], we first find a Grobner basis G 
for V for the lexicographic order yo > xq > x > y and reduce it eventually to a unique reduced basis G r of 
nine (or eleven, depending on the case) polynomials. Finally, we eliminate from that basis all polynomials 
that contain variables yo and xq, and keep the single polynomial g G R[x, y] that gives the parallel lines and 
generates the second elimination ideal fa = (g)- We will also proceed to find all singular points of V(g) where 
(Vg)(p) = 0. Thus, let's consider three cases when < r < r crit , r = r crit , and r > r crit . In the following 
section we will consider a general case for arbitrary values of a, b (a > b > 0) and r (r > 0). 
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Example 7 (r < r cr i t ) Let a = 3, b = |, and r = ^. Then, 

v = v (yo -x"o-l, 4y 2 - 8yy + 4y 2 + 4a; 2 - 8a;a;o + 4a; 2 , - 1, xy + yx - 2y Q x Q ) 
We find that the reduced Grobner basis G r contains nine polynomials including 

g = 225 + 64/ + 288a; 6 + 480a; 2 ?/ 4 - 800yV - 480y 6 - 128y 4 x 4 

- 1180y 2 + 516a; 2 - 1320a; 2 y 2 + 532a; 4 + 1236y 4 + 64a; 8 (52) 

whose gradient is 

Vg = a;(216x 4 + 120y 4 - 400x 2 y 2 - 64x 2 y 4 + 129 - 330y 2 + 266a; 2 + 64a; 6 )i 

(53) 

- y(-64y 6 - 240x 2 y 2 + 200a; 4 + 360y 4 + 64y 2 a; 4 + 295 + 330a; 2 - 618y 2 )j . 
There are exactly two solutions to the equation (V<?)(p) — on V(g) : 

We will refer to these points as virtual singular points for the following reason: V(y) is a disjoint union 
of {Si, S2} and a non-singular smooth subvariety II where the gradient V<? does not vanish. The subvariety 
II constitutes the set of two parallel lines that we have been seeking. Note that the foci F\ = (0, y/2) and 
F2 = (0, — v2) are different from the singular points. We show all in Figure^ 

We have verified by direct computation that as the value of the offset r increases towards some critical 
value r cr it, the virtual singular points approach the smooth component of the variety V(^), that is, they 
approach the parallel lines II. 7 As the next example shows, when r = r cr n, then the virtual singular points 
are actually located on the smooth component II of the variety V(y). 

Example 8 (r = r cr i t ) Let a — |, b — 1 , and r — | . Then, 

V = V(4y 2 - 9a; 2 , - 9, 9y 2 - I8yy + 9y 2 + 9x 2 - 18xx a + 9x 2 ~ 4, 4a;y - 13a; y + 9a;oy) 

We find that the reduced Grobner basis G r contains eleven polynomials including 

g = 9447840a;V - 4199040y 6 x 2 + 53800200a; 6 + 1679616y 8 + 124857369x 4 

+ 127471968?/ - 159339960x 2 y 2 - 24820992y 6 + 46539360yV - 105471720y 2 a; 4 

- 4933872y 4 a; 4 + 120670225 + 8503056a; 8 - 250879824y 2 + 156799890a; 2 (54) 

whose gradient is 

V.9 = a;(1574640y 2 a; 4 - 233280y 6 + 8966700a; 4 + 13873041a; 2 - 8852220y 2 + 2585520y 4 

- 11719080a; 2 y 2 - 548208yV + 1889568a; 6 + 8711105)1 

(55) 

+ y(131220a; 6 - 174960y 4 a; 2 + 93312y 6 + 3540888y 2 - 2213055a; 2 - 1034208y 4 

+ 1292760a; 2 y 2 - 1464885a; 4 - 137052y 2 a; 4 - 3484442)j . 

7 In the next subsection we will derive a formula for the critical value r cr i t of r for a general hyperbola. 
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Figure 7: Hyperbola with parallel lines <?, foci F\ 1 F2 and two singular points Si, S3 for a = 1,6 = 1 and 

r = i 
' 2 



There are exactly two solutions to the equation (X7g)(p) — on *V(g) : 

St = (o,^f\ and s 2 = (o,~ 

These points are the singular points of the variety "V(g) : At every other point of this variety the gradient 
V<7 does not vanish. We will collect these points as before into a non-singular subvariety II. The variety 
II constitutes the set of two parallel lines that we have been seeking. Note that the foci F\ = (0, and 
F2 = (0, — ^p-) are still different from the singular points. We show all in Figure^ 

As we can see from the last figure, when r — r cr u — | for the given values of a and 6, the singular points 
are located on the continuous component of II. 

Probably the most interesting parallel lines we again obtain when r > r cr i t as we show in the next 
example. There, we will find six singular points as each of the two singular points found in Example |S] splits 
into three points. 

Example 9 (r > r cr i t ) Let a = |, b = 1, and r = |. Then, 

V = V(4y2 - - 9, 9y 2 - 18 yy() + 9yg + 9.t 2 - 18xx + 9x% - 16, 4xy - 13x Q y + 9x y) 
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Figure 8: Hyperbola with parallel lines g, foci Fi,F 2 and singular points Si, S2 for a = §, b = 1 and r = | 



We /md tftat i/ie reduced Grobner basis G r again contains nine polynomials including 

g = -1627128a; 6 + 8503056a; 8 + 30525625 + 83351034a; 2 - 1109487600y 2 + 264529800a; 2 y 2 

+ 51521913a; 4 + 381560544y 4 + 349920y 4 a; 2 - 188052840yV - 43856640y 6 - 4933872yV 



Vg = x(-271188x 4 + 1889568a; 6 + 4630613 + 14696100y 2 + 5724657a; 2 + 19440y 4 

- 20894760xV - 548208y 4 x 2 + 1574640j/V - 233280y 6 )i 

(57) 

+ y(-15409550 + 3674025x 2 + 10598904y 2 + 9720a; 2 y 2 - 2611845a; 4 - 1827360y 4 

- 137052j/ 2 a; 4 + 131220a; 6 - 174960j/ 4 a; 2 + 93312y 6 )j . 

There are exactly six solutions to the equation (Vg)(p) = on V(g), or, in another words, six singular points 
of the variety: 



+ 9447840a;y - 4199040yV + 1679616y 8 (56) 



whose gradient is 
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'2^39- 78-21? + 39 • 2§ ^12805 + 11232 ■ 2s + 12636 • 2i 



Sa=\- 



13 78 
2^39- 78-23 +39-2§ a/12805 + 11232 • 2s + 12636 • 2r 



13 



78 



Sn = 



2 a/39 - 78 • 2? + 39 • 2§ \/l2805 + 11232 • 2s + 12636 • 2i 



13 



78 



5 B = - 



2^39- 78-23 +39-2i a/12805 + 11232 • 2? + 12636 • 2? 



13 



78 



We display these six points in Figure^ 




Figure 9: Hyperbola with parallel lines g, foci Fi, F2 and six singular points Si, i — 1, ... , 6, for a = |, b = 1 
and r = I 



4.2 Parallel lines for arbitrary parameters a, 6, r 

In this section we will find a general polynomial g that will generate the second elimination ideal 1% = 
I n M\x, y, a, b, r] where I = /a, /a) and the generating polynomials are defined in equations l|i9"|l. 
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and 151|) . We will also find a formula for the critical value r cr n of r that determines whether V(i2) has two 
or six singular points, and, we will find general formulas for the coordinates of the said singular points for 
arbitrary values of the parameters a, b and r. 

We begin with a computation of a Grobner basis for the ideal / for the lexicographical order yo > xq > 
x>y>a>b>r. 

I = (6 2 y 2 - a 2 x% - a 2 b 2 ,y 2 - 2yy + y 2 + x 2 - 2xx + x\ - r 2 ,b 2 y a x - b 2 y a xa + a 2 x y - a 2 x y ) 

We found that a reduced Grobner basis for / consists of fifteen homogeneous polynomials of degrees: 16 
(two polynomials), 14 (one polynomial), 12 (one polynomial), 11 (two polynomials), 10 (two polynomials), 
9 (one polynomial), 8 (one polynomial), 7 (two polynomials), 5 (one polynomial), 4 (one polynomial), and 
2 (one polynomial). We will not display these polynomials as they are similar to the polynomials for the 
ellipse shown in Appendix B. 

Thus, in this general case we observe that I2 = (g) where g = g\/{a 2 b 2 ) is the following homogeneous 
polynomial in M.[x, y, a, 6, r] of degree twelve: 

g = 6 a 4 x 4 r 4 - 2 b 2 x 4 a 2 r 2 y 2 + 4 b 6 x 2 a 4 + 2 a 8 b 2 x 2 + 2 a 6 x 6 + 6 y 4 a 2 b 6 - 2 a 4 b 2 x 4 y 2 
+ 4a 4 x 6 b 2 -6a e x 2 b 2 y 2 - 4 a e x 2 b 2 r 2 - 10 a 4 x 2 b 4 y 2 - 8 a 4 x 2 b 4 r 2 
+ 8a 4 x 2 b 2 r 4 -6 b 4 r 2 a 2 x 4 + 4 b 4 r 4 a 2 x 2 - 10 a 4 x 4 b 2 r 2 - 6 b 4 r 2 a 2 x 2 y 2 

1 a 4 2 7.2 2 2 £.7.4 4 2 2i<-> 7,4 2 4 2 ei.6 222iinj.4 422 

+ 6 a x b r y — b x a y + 2 x y a — 6 x r a + 10 y a r 

o 7,4 2 4 2 1 1,6 222 o 1,4 242,^44 i2 2 o „2 2 t8 , c i4 4 4 

- 8 b y a r — 4b y r a — 8 b y r a + by a b x — 2a y o +00 y a 

- 4 b 4 y 6 a 2 - 6 y 2 a 4 b 6 + 6 a 4 x 4 6 4 - 4 y 2 a 6 6 4 - 6 y 2 a 6 r 2 b 2 + 6 y 4 a 4 r 2 6 2 

- 4 y 2 a 4 r 4 rb 2 + 2 a 6 6 6 + a 4 6 8 + 6 a 6 6 4 x 2 + 6 a 6 x 4 b 2 + b s r 4 + a s b 4 - 2 a 2 b s r 2 

- 2 a 2 b 6 r 4 + 2 b 6 x 2 r 4 - 4 b 6 x 2 y 2 r 2 + 6 b 2 x 4 a 2 r 4 + 4 a 6 x 2 r 2 y 2 - 2 fb 2 x 2 a 2 y 6 

+ 2b 4 x 2 y 6 -2b 2 x 2 y 4 a 2 r 2 + 6 5 4 x 2 y 2 r 4 + 10 6 2 x 2 y 2 a 2 r 4 + 2 6 6 x 2 y 4 (58) 

- 6 b 4 x 2 y 4 r 2 +a 4 x 8 + 2 b 6 r 6 + b 4 r 8 + 6 a 6 x 2 r 4 - 6 a 6 x 4 r 2 - 2 a 8 x 2 r 2 

- 6 a 4 b 4 r 4 - 2 a 4 b 6 r 2 - 2 a 6 y 2 r 4 - 2 a 6 b 2 r 4 - 6 a 2 r 6 6 2 y 2 - 2 a 4 r 6 y 2 

- 6 a 4 x 4 r 2 y 2 + 6 a 2 y 4 b 2 r 4 + 2 a 8 b 2 r 2 + a 4 y 4 r 4 - 2 a 2 y 6 b 2 r 2 - 4 a 4 r 6 x 2 
+ a 8 x 4 + 6 a 4 x 2 y 2 r 4 - 4 a 4 x 6 r 2 + a 8 r 4 + 2 a 4 x 6 y 2 - 2 a 6 x 4 y 2 - 2 a 6 r 6 
-2a 4 r e b 2 + 2a 6 r 2 b 4 + a 4 y 4 x 4 + 2 a 2 r 8 b 2 - 2 a 4 y 4 x 2 r 2 + 2 a 2 b 4 r 6 + a 4 r 8 

+ 6 y 4 b 4 r 4 + 6 y 4 b 6 r 2 - 4 y 6 6 4 r 2 - 4 y 2 6 4 r 6 - 6 y 2 6 6 r 4 + y 4 b 8 + y 8 b 4 - 2 y 6 6 6 

ni.2 622 o 7,4 26 ot.2 262,7.4 44 it.2 442,7.4 44 

- 2b x a y —2bxr —obxr a +b x y — 4b x y a +b x r 

o 7,2 6 2 2 o 7,4 422 r<2 7,8 2 a 2 2 2 iS 

- 2b x a r —2 b x r y — 2 y o r — 6 y x a o 

We want to find now those points on the variety V = V(g) where the gradient Vy is zero. Thus, V can 
be thought of as union of points where the gradient is zero and those points where the gradient is not zero. 
Connected components of V where the gradient is not zero constitute the 1-dimensional manifold. We are 
interested in a discrete finite set S of points where the gradient is zero. These points will give the singular 
points of the variety V. We will show next that the critical value of r is r cr ;t = — ■ Therefore, we compute 
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= Vg(x, y, a, b, r) =< h\, hi, /13, h 4 , /i 5 > where 

h 1 =x(2b 6 y 2 r 2 -a s b 2 + b i r 6 -x 2 b 4 r 4 + 2x 2 b 4 r 2 y 2 - 3b 4 y 2 r 4 + 3 6 4 y 4 r 2 - b 6 y 4 

- 3 a 6 b 4 - a 8 x 2 - 2 a 4 b e - 3 a 4 y 4 b 2 + 3 y 2 a 6 o 2 + 3 y 2 o 6 a 2 + 5 o 4 y 2 a 4 

- 6 4 y 4 a 2 — 6 a 4 x 2 o 4 + r 2 a 8 - o 4 x 2 y 4 + 2 a 6 x 2 y 2 + 6 o 4 x 2 a 2 y 2 + 2 a 4 o 2 x 2 y 2 

- 6 4 y 6 - b 6 r 4 + 10 a 4 x 2 o 2 r 2 - 6 a 6 x 2 b 2 - 6 a 4 x 4 & 2 - 3 a 6 x 4 + 3 o 6 r 2 a 2 

- 2 o 4 r 4 a 2 + 2 r 6 a 4 + 3 r 6 o 2 a 2 - 3 r 4 a 6 + 4 r 2 a 4 6 4 + 2 b 2 r 2 a 2 x 2 y 2 

+ 6 a 4 x 2 r 2 y 2 + 6 a 6 x 2 r 2 + 3 b 2 r 2 a 2 x 4 - 6 a 4 x 2 r 4 + o 2 y 4 a 2 r 2 - 5 o 2 y 2 r 4 a 2 

o244 ci.2 4 2 2ioi.2 4 2 2iyii.2 2 4 2 442 o442 

- 3y a r —do r a x + 3 b x a y + A b x y a — y a x — 3 a x y 
+ 6 a 4 x 4 r 2 - 2 y 2 a 6 r 2 + y 4 a 4 r 2 + 6 2 y 6 a 2 + 3 o 4 y 2 a 2 r 2 - 3 & 2 y 2 a 4 r 2 

+ 2o 2 a 6 r 2 +6o 4 x 2 a 2 r 2 -4o 2 a 4 r 4 -2a 4 x 6 ) 

h 2 = y (-6 6 6 y 2 r 2 + 2 6 4 r 6 - 3 x 2 o 4 r 4 + 6 x 2 6 4 r 2 y 2 - 6 6 4 y 2 r 4 + 6 o 4 y 4 r 2 + r 2 b 8 
+ b 4 x 4 r 2 + 2 x 2 6 6 r 2 + 3 6 6 y 4 + 2 a 6 o 4 + 3 a 4 6 6 - 6 y 2 6 6 a 2 - 6 6 4 y 2 a 4 
+ 6 b 4 y 4 a 2 + o 2 x 6 a 2 + 3 o 6 x 2 a 2 + 3 o 4 a 2 x 4 + 5 a 4 x 2 6 4 + a 2 6 8 - 2 y 2 x 2 6 6 

- 3 6 4 x 2 y 4 - o 4 x 4 y 2 - 2 6 4 x 2 a 2 y 2 - 6 a 4 o 2 x 2 y 2 - 2 o 4 y 6 + 3 o 6 r 4 

- 3 a 4 x 2 b 2 r 2 - y 2 6 8 + 3 a 6 x 2 6 2 + a 4 x 4 6 2 + a 6 x 4 + 2 b 6 r 2 a 2 + 4 6 4 r 4 a 2 + r 6 a 4 

io6l2 2, 461,124 1.4 10 1,2 2222104222 o 6 2 2 

+ 3 r o a + r a +4r a o +26 r a x y +2a x r y — 2a x r 

.2 2 2 4 o424iol.2 422 c 1.2 2 4 2 244 r i2 4 2 2 

+ o r a x — 3 a x r + 3 o y a r — 6 o y r a — y a r — o b r a x 
+ 4b 2 x 4 a 2 y 2 + 3 b 2 x 2 y 4 a 2 - a 4 x 4 y 2 + 3 a 4 x 4 r 2 - 10 b 4 y 2 a 2 r 2 - 6 6 2 y 2 a 4 r 2 
+ 36 2 a 6 r 2 + 36 4 x 2 a 2 r 2 + 2 6 2 a 4 r 4 - a 4 x 6 ) 

h 3 = 2b 2 x 4 y 4 + 6 2 x 6 y 2 + 6 2 x 2 y 6 + 2 6 6 y 2 r 2 - 6 4 r 6 - 2 x 2 6 4 r 4 + 3 x 2 6 4 r 2 y 2 

+ 4 6 4 y 2 r 4 - 5 6 4 y 4 r 2 + r 2 6 8 + 3 b 4 x 4 r 2 + 3 x 2 6 6 r 2 - 3 6 6 y 4 - 2 a 6 6 4 - a 2 x 8 
+ b 2 x 4 r 2 y 2 + o 2 x 2 y 4 r 2 - 5 b 2 x 2 y 2 r 4 - a 2 y 4 x 4 - r 8 b 2 + 2 a 2 r 6 y 2 - a 2 y 4 r 4 
+ 4 a 2 r 6 x 2 - 6 a 2 x 4 r 4 - 3 o 2 x 4 r 4 + 6 a 2 x 4 r 2 y 2 - 3 y 4 6 2 r 4 + 3 r 6 6 2 y 2 + y 6 6 2 r 2 

c2224 l r ) 2422 l0 J,2 26 1 j.2 62 l/ |262 r>262 82 

- 6 a x y r +2a y x r +3o x r + 6 x r +4a x r — 2a x y — r a 

- 3 a 4 6 6 + 6 y 2 o 6 a 2 + 6 6 4 y 2 a 4 - 6 6 4 y 4 a 2 - 4 o 2 x 6 a 2 - 4 6 6 x 2 a 2 - 6 o 4 a 2 x 4 

- 9 a 4 x 2 6 4 - a 2 6 8 + 3 y 2 x 2 b 6 - b 4 x 2 y 4 + 3 o 4 x 4 y 2 + 10 b 4 x 2 a 2 y 2 

+ 9a 4 b 2 x 2 y 2 + 2 b 4 y 6 + 6 6 r 4 + 6 a 4 x 2 b 2 r 2 + y 2 6 8 - 4 a 6 x 2 6 2 - 9 a 4 x 4 o 2 
-2a 6 x 4 + 2o 6 r 2 a 2 + 66 4 r 4 a 2 + 3r 6 a 4 + 2r 6 o 2 a 2 - 2 r 4 a 6 - 3 r 2 a 4 6 4 

- 6 b 2 r 2 a 2 x 2 y 2 - 6 a 4 x 2 r 2 y 2 + 4 a 6 x 2 r 2 + 10 6 2 r 2 a 2 x 4 - 9 a 4 x 2 r 4 

- 6 b 2 y 4 a 2 r 2 + 4 6 2 y 2 r 4 a 2 + 3 y 2 a 4 r 4 - 8 6 2 r 4 a 2 x 2 + 2 b 2 x 4 a 2 y 2 

- 6 b 2 x 2 y 4 a 2 + 3 a 4 x 4 y 2 + 9 a 4 x 4 r 2 + 8 6 4 y 2 a 2 r 2 + 9 o 2 y 2 a 4 r 2 - 4 6 2 a 6 r 2 
+ 86 4 x 2 a 2 r 2 + 36 2 a 4 r 4 -3a 4 x 6 
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K = b 2 x 4 y 4 + 2 b 2 x 2 y^-A 6 6 y 2 r 2 + a 8 b 2 + 3 6 4 r 6 + 3 x 2 b 4 r 4 - 6 x 2 b 4 r 2 y 2 

-9b 4 y 2 r 4 + 9b 4 y 4 r 2 + 2b 6 y 4 + 3a 6 b 4 + b 2 y 8 + a 8 x 2 ~2b 2 x 4 r 2 y 2 

— 66 2 x 2 y 4 r 2 +6b 2 x 2 y 2 r 4 - 2 a 2 y 4 x 4 + r s b 2 -3a 2 r 6 y 2 + 3 a 2 y 4 r 4 - 3 a 2 r 6 x 2 

, o 2 4 4 , j,2 4 4 2422,c4i2 4 ,i 6 7,2 2 a 6 7,2 „2 , c 2 2 2 „4 

+ 3(1 fr I o x r — a x r y + 6 y or — A r by — 4 y b r +5 a x y r 

2422 o 1.2 26 262 262,82 226 2 6 2 , n 4 7,6 

— a y x r — 2 o x r — a x r — a x y + r a — x a y — a y r + 2 a 6 
+ 6 a 4 y 4 6 2 - 4 y 2 a 6 6 2 - 4 y 2 b 6 a 2 - 9 & 4 y 2 a 4 + 9 6 4 y 4 a 2 + 6 a 4 x 2 b 4 + r 2 a 8 

+ 3 b 4 x 2 y 4 - 3 a 6 x 2 y 2 - 9 b 4 x 2 a 2 y 2 - 10 a 4 6 2 x 2 y 2 - 3 6 4 y 6 + 2 6 6 r 4 (62) 

— 8 a 4 x 2 6 2 r 2 + 6 a 6 x 2 fe 2 + 6 a 4 x 4 6 2 + 3 a 6 x 4 - 4 & 6 r 2 a 2 - 3 6 4 r 4 a 2 - r 6 a 4 
+ 2 r 6 6 2 a 2 - r 4 a 6 - 3 r 2 a 4 6 4 - 6 6 2 r 2 a 2 x 2 y 2 + 3 a 4 x 2 r 2 y 2 - 2 a 6 x 2 r 2 

e i2 2 2 4 , j 4 „2 „4 , 1 r> 1.2 „ 4 2 „2 oj.2 242 r,244,^ 7.2 422 

— 6o r a x +4axr + 10 b y a r —sbyra — 2y a r +46 r a x 

rJ,2 422ioJ,2 242io442 442 r442 o262 

— 6 o x a y + 2 b x y a + 3 y a x — a x y — 5 a x r — 3 y a r 

,o442 ,< 1,2 62 a lA 2 2 2 oj.2 2 4 2,0 .2 6 2 n 1.4 222 

+ 3 y a r — 4 b y a — 6 o y a r — 8 b y a r +2 b a r — 9 o x a r 
-66 2 a 4 r 4 + 2a 4 x 6 



(63) 



h 5 = -6b 6 y 2 r 2 +a s b 2 + 2b 4 r 6 -3x 2 b 4 r 4 + 6x 2 b 4 r 2 y 2 - 6 b 4 y 2 r 4 + 6 b 4 y 4 r 2 + r 2 6 8 
+ b 4 x 4 r 2 + 2 x 2 b 6 r 2 + 3 6 6 y 4 + a 6 b 4 - a 8 x 2 - a 4 b 6 + 3 a 4 y 4 6 2 - 3 y 2 a 6 6 2 

- 2 y 2 fe 6 a 2 - 4 6 4 y 2 a 4 + 5 b 4 y 4 a 2 - fe 2 x 6 a 2 - 3 fe 6 x 2 a 2 - 3 6 4 a 2 x 4 - 4 a 4 x 2 6 4 

, 28 2 j.8 o22i6 o 1.4 24 !.4 42,o622 o 1,4 2 2 2 

+ r a — a b — 2 y x o — 3 o x y — b x y +2 a x y — 3 o x a y 

,o4t2 22 oi.4 6,o;.6 4io42i2 2 2 18 o62j,2 r44 7,2 

+ 3a x y — 2 o y + So r + 8a x b r — y b —2a x b — 5 a x b 

- 3a 6 x 4 - 2& 6 r 2 a 2 + 3 6 4 r 4 a 2 + 2r 6 a 4 + 4r 6 fe 2 a 2 - 3 r 4 a 6 - 6 r 2 a 4 b 4 
+ 10 6 2 r 2 a 2 x 2 y 2 + 6 a 4 x 2 r 2 y 2 + 6 a 6 x 2 r 2 + 6 b 2 r 2 a 2 x 4 - 6 a 4 x 2 r 4 

,/3 (.2 422 n 7.2 2 4 2 o244 nl.2 422 t.2 422 j.2 242 

+ 6 o y a r — 9 o y r a — 3 y o r — 9 b r a x — b x a y — b x y a 

442 o442,c442 r,262, 442 t.2 62 o 1,4 2 2 2 

— y a x — 3 a x y + 6 a x r — 2 y a r + y a r — b y a — 8 b y a r 

a 7.2 2 4 2 oj,2 6 2 ly( i,4 2 2 2 o!,2 44 r.46 

— Ab y a r —2b a r +46 x a r — 36 a r — 2a x 

Continuing our approach from Section [21 in order to find the solutions to the vector equation Vg = 0, we 
need to solve a system of six polynomial equations hi = 0,i = 1, . . . , 5 and g = in the variables x, y, a, 6, r 
under the assumption that a > b > r > 0. Let Is — (Vg,g) — (hi,h2,h$,hi,h5,g) and let F be a list 
containing the six polynomial equations: 

F = [hi = 0, h 2 = 0,h 3 = 0, hi =0,h 5 = 0, g = 0] (64) 

In order to solve system (|64() . we can again either employ the Grobner basis approach similar to the one 
shown in Appendix C in case of the ellipse, or simply use solve command in Maple. Thus we find that the 
affine variety V(7s) consists of the following eight singular points of which six are real and two are complex: 



Si= L W + »Ur*+n), ga= fo,.vg±S£±gn, (65) 
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S3 = Mg±gEg±gU, a,-f-^±gEg±gI,o|. (66) 



£ y / (r 6) 3 (a 2 + 6 2 )e 



^5 = « - 7 „ , — -r , / tttt^t - ^ - (67) 



^6= -J — - , —7 , n2 ,- V , , , (68) 



(a 2 + b 2 )(r b)i ' 


(r b) i (a 2 H 


-6 2 ) y 
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S7=\ o . ,„w „a > V , 7TT7 =— , (69) 



(a 2 + 6 2 )(r6)§' (r 6) 3 ( a 2 + 6 2 ) 



(70) 



where 



S = 3 6 2 J r 2 - r 2 a 2 (r b) § - 3 a* (r 6) (1/3) 6 3 r + 6 4 (r 6) i , (71) 
£ = a 4 (r6)i +6 2 r 2 (r£>)s +3 6a3 r (r&)s + 3cJ 6r. (72) 
Notice that <5 becomes zero in three cases: 

1. When r = - this is not possible since we expect r > 0, 

,2 

2. When r = — — - this is not possible as r > 0, 

3. When r = ^. 

a 

Only this last value is possible thus we set r C rii = — • When r = r cr u 5 then under the assumptions that 
a > b > r > 0, we essentially obtain two real singular points: 

Si = S 7 = S s = (o, = (o, and S 2 = S 5 = S 6 = U ~^^) = (o, -£) , (73) 

while the complex points S3 and S4 become: 

S 3 = + & Va 2 -6 2 i i 0^ and 5 4 = (' "'^V a 2 - 5 2 »,oj ■ (74) 

,2 

For a hyperbola with a > b, the ratio i = — is called a semi-latus rectum. It is the distance between a 
focus and the hyperbola itself measured along a line perpendicular to the major axis. The value of / is also 
the reciprocal of the maximum curvature K max of the hyperbola, or, it is equal to the radius of the smallest 
osculating circle. \TQ 

To summarize, the variety 'V(Is) contains the following real singular points: 



27 



1. Two virtual points S\ and S2 shown in (|65f) when < r < r cr n. In particular, when a = 1, b = 1, and 
r = i, upon substituting these values into l|58[) we obtain the parallel lines (|52(l derived in Example [7] 
Furthermore, when these values are substituted into I|tj5[) . then we obtain the two virtual singular 
points Si and S2 found in the same example. The remaining six singular points are complex. 

2. Two points Si and £2 shown in (|65[1 when r — r cr u. In particular, we can recover results of ExampleEl 
when we substitute values a = §, b = 1, and r — | into the general equation 1)58(1 and the singular 
points J73J). To be precise, three of the singular points give Si, another three give S2 whereas the 
remaining two gave two distinct complex points. 

3. Six points Si, S2, S5, Sq, S7, S$ shown in jBSjl , (jBTJl, JBSJ, and ifTUjl when r > r cri t. In particular we 
can recover in the same manner six points and the equation for the parallel lines found in Example 
when we set a = |, b = 1, and r = |. These six singular points are displayed in Figure El The 
remaining two points are complex and distinct. 

5 Applications and Conclusions 

There are several immediate applications for the results obtained in the previous sections as well as for the 
general method of solving systems of polynomial and rational equations using the Grobner basis technique. 
In particular, we refer to the explicit form of the polynomial g that defines the parallel lines to the conies 
given in (|T^|) . i|32[l. and While finding g amounts to computing the reduced Grobner bases for the 

appropriate ideals and then applying the elimination theory - both being standard techniques of algebraic 
geometry 0], 0, [E], [E], H2| ~ these techniques are not used and probably not even known in the area of 
engineering applications where numeric, CAD generated curves, dominate. However, it does not mean that 
the fact that the critical value r cr n of r equals the reciprocal of the maximum curvature of the curve - as 
derived above for the non-generate conies - is not known in engineering: In fact, this fact is very well known 
in the theory of mechanisms when designing cams with reciprocating roller follower [111 Chapter 8], |13l 
Section 5.10]. 

In the theory of machines, the "inner" parallel line (on the concave side of the curve) is referred to as 
the cam profile, while the path followed by a center of a roller - in our case it was always given above as 
the polynomial fi - is referred to as the pitch curve |13|. The pitch curve is usually a smooth curve 
defined piece-wise by the pre-determined motion of the follower that could be harmonic, or determined by a 
polynomial of degree eight (the so called displacement function) . This results in more complicated - than in 
our examples above - cam contours that may have portions which are concave, convex or flat jlll page 404]. 
The problem of the inner curve being "non smooth", or, in our language, having singularity points on the 
smooth portion of the variety V(g), is referred to as undercutting - when the radius r of the roller follower is 
greater than r cr i t . When r = r cr n , that is when the radius of the roller follower equals the minimum positive 
(convex) radius of the pitch curve, a cutter cutting out the cam will create a cusp point on the cam surface. 
In this latter case, while the cam contour is a continuous curve, it is not smooth at one or more points that 
results in a cam not running well at high speeds. In the first case when r > r cr i tl the cutter cutting the 
cam out (...) undercuts or removes material needed for cam contours in different locations and also creates 
a sharp point or cusp on the cam surface. This cam no longer has the same displacement function you so 
carefully designed jlll page 404]. Such cams are not acceptable as the radius of the roller follower is too 
large. As a rule of thumb, the absolute value of the minimum radius of curvature of the cam pitch curve is 
preferably at least 2 or 3 time as large as the radius of the roller follower r. These situations that correspond 
to our cases discussed in examples El and 01 for the parabola, El and El for the ellipse, and |H1 and El for the 
hyperbola, are depicted very well in ^2 Figure 8-45] and in ^3 Figure 5-29]. 
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What differs in the approaches to the problem of parallel lines presented above and designing cam profiles 
is that (i) in the former approach we obtain explicit formulas for the parallel lines given by the polynomial g 
rather than obtain a numeric representation for it using some CAD software, (ii) in reality cam profiles are 
more complicated as they are not determined by a single curve but are determined piecewise and, preferably, 
smooth-wise. Our approach allows for an analytic, hence exact, representation of the parallel lines, and it 
therefore affords an exact differential calculus to be applied to g. There is no reason why in principle this 
approach could not be extended to piece-wise defined pitch curves. In particular, it allows for the exact 
computation of the coordinates of the singular points. 

In the following, we briefly show one application of the result presented in Example^] Finally, we remark 
that these methods extend to non conies such as a Bezier cubic 

5.1 Parallel lines to an ellipse in the finite element method 

Let us consider an ellipse determined by a polynomial fx : 

fi =4^ + 162^-64 (75) 

for the values of a = 4 and 6 = 2. Thus, c = 2-\/3 and r cr n = — = 1. We will generate a finite element mesh 
by choosing these values for the offset r = 0.2, 0.4, 0.6. These numbers determine thicknesses of layers, or, 
the distances of the parallel lines to the ellipse. Number of rows in the mesh matrix is twice the number of 
values of r values plus 1. 

Remark 1 The reason we are interested in generating several smooth parallel lines to the ellipse is that we 
want to apply the finite element analysis to composite materials that are built precisely from smooth layers of 
various materials. As long as the parallel lines are smooth, and this is guaranteed in the range < r < r cr n, 
the layers are smooth and touch at each point. This guarantees, in turn, an optimal performance of the 
composite. 

For each of the chosen values of r, we can now compute the corresponding polynomial g by substituting 
values of a,b,r into Thus, we get: 

,9i = -0.1083369141 • 10 8 y 2 - 0.1281707488 • 10 8 x 2 + 0.4086271181 • 10 7 y 2 x 2 

+ 33792 y 4 x 4 + 10240 x 2 y 6 + 40960 x 6 y 2 - 495206.40 x 4 y 2 - 369868.80 x 2 y 4 

+ 259850.24 x 6 - 57180.16 y 6 + 0.1185870643 • 10 7 y 4 - 564800.7168 x 4 + 1024 y 8 

+ 16384 a; 8 + 0.3681277971 • 10 8 when r = 0.2, (76) 

.92 = -0.1031174671 • 10 8 y 2 - 0.1349849501 • 10 8 x 2 + 0.4156607693 • 10 7 y 2 x 2 

+ 33792 y 4 x 4 + 10240 x 2 y 6 + 40960 x 6 y 2 - 506265.60 x 4 y 2 - 373555.20 x 2 y 4 

+ 252968.96 a; 6 - 56688.64 if + 0.1155347251 • 10 V - 684903.6288 a; 4 + 1024 y 8 

+ 16384 a- 8 + 0.3409692932 • 10 8 when r = 0.4, (77) 

.93 = -0.9465872249 • 10 7 y 2 - 0.1456495775 • 10 s x 2 + 0.4277898445 • 10 7 y 2 x 2 

+ 33792 y 4 x 4 + 10240 x 2 y 6 + 40960 x 6 y 2 - 524697.60 x 4 y 2 - 379699.20 x 2 y 4 

+ 241500.16 x 6 - 55869.44 y 6 + 0.1104343859 • 10 7 y 4 - 880291.0208 a; 4 

+ 0.2986886575 • 10 8 + 1024y 8 + 16384 x 8 when r = 0.6. (78) 

Thus, so far we have rows in a mesh matrix. In order to create columns, we select some y values that 
determine the y-coordinates of points on the ellipse through which we draw perpendicular lines to the 
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ellipse. The intersections of these lines with the ellipse itself as well as with the parallel lines are the mesh 
nodes. The number of these y- values gives the number of mesh columns minus 1. The last column is obtained 
by including intersection points between gx, g2, 53, and f± with the y-axis. For the sake of this example, we 
select these values: y = 3.75, 3, 2, 1,0,— 1,-2,— 3,— 3.75 (the y- values must satisfy —a < y < a). We display 
the ellipse, the parallel lines, and the nodes in Figure ITU1 
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Figure 10: Ellipse with parallel lines g\, g%, 53 and nodes for a = 4, b = 2 and r = 0.2, 0.4, 0.6 

It is easy now to produce a 7 x 20 mesh matrix with 140 nodes and to create a list of 120 4-node elements 
and a list of 30 9-node elements. Once the elements have been created in Maple, output has been sent to 
a finite element software for analysis of natural frequencies. It has been found, that there is a significant 
sensitivity in the higher mode of natural frequencies. For instance, when the mesh generated by exact parallel 
lines has been used, the higher mode of natural frequencies are much higher than the ones obtained from a 
mesh generated based on approximately parallel lines - several confocal ellipses. 

A Appendix: Reduced Grobner basis for the ideal 

I = (fh /2, h) C R[y 0: x 0: x, y, r,p] 

For the ideal 

I = (Apy Q - xl,y 2 - 2yy + y\ + x 2 - 2xxq + x 2 - r 2 ,2xp- 2x p + x y - x y Q } C R[y ,x ,x,y,r,p] 
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using lex order with y > x a > x > y > r > p, a reduced Grobner basis for I is given by the following 
fourteen polynomials: 8 

> gl :=f actor (reducedG [1] ) ; 

gl := p(-2pr 2 yx 2 + 8pr 2 y 3 + 8p 2 r 2 y 2 - 32 yp 3 r 2 + 16 p 4 r 2 - 16 y 4 p 2 + 32y 3 p 3 

- 16p 4 y 2 + 3r 2 x 4 + 8p 2 r 4 + 20p 2 r 2 x 2 - y 2 x 4 + lOypx 4 - x e - x A p 2 + 8py 3 x 2 

- 32 x 2 y 2 p 2 + 8x 2 yp 3 - 3 r 4 x 2 + 2 r 2 x 2 y 2 + r 6 - r 4 y 2 - 8pr 4 y) 

> g2 : =f actor (reducedG [2] ) ; 

g2 := p(-lUx 3 p 2 y + 232 xp 3 y 2 + 40 xp 4 y- 316 xp 3 r 2 - 8xy 3 p 2 - 2x 3 p 3 + 16 x p 5 

- 16 py 2 xr 2 - \2py 2 xO r 2 + 54xyp 2 r 2 - 120 yp 2 xO r 2 - 32 yp 4 xO 

+ 104 p 3 xO r 2 - 14py 2 x 3 + 32p 2 y 3 xO + 32xpr 4 - 46x 3 pr 2 - 16p 5 xO 

+ 14 x 5 p + 27 r 4 xO p + 2 y x 5 + 2 y 3 x 3 - 2 y 3 x r 2 - 4 y 3 xO r 2 + 2 y x r 4 - 4 y x 3 r 2 

- 24 py 4 x + 16 py 4 xO) 

> g3 : =f actor (reducedG [3] ) ; 

g3 := p(4x 2 p 3 + 56 p 3 r 2 - 4xp 3 xO - 56 p 3 y 2 + 48 y 3 p 2 + 42x 2 p 2 y - 12 xp 2 yx0 

- 48 yp 2 r 2 — 4x 4 p + 8 y 4 p + 14 pr 4 — 10 x 2 pr 2 — 12 xp y 2 xO — 8x 2 py 2 

+ 27 xp xO r 2 - 22 y 2 p r 2 + 2 r 4 y + 2 y 3 x 2 + 2 y x 4 - 4 r 2 y x 2 - 4 x y 3 xO - 2 r 2 y 3 ) 

> g4 : =f actor (reducedG [4] ) ; 

g4 := p(112xp 4 - 112x0p 4 - 120yxp 3 + Y16yx0p 3 - 124xp 2 r 2 - 80y 2 x0p 2 

- 16 xOp 2 r 2 + 160 x p 2 y 2 - 96x 2 p 2 xO + 82 x 3 p 2 - 28 x 3 y p + 4 y xOpr 2 

- 24 pxy 3 + 16 py 3 xO + 22 xypr 2 - 4r 2 y 2 xO + 2y 2 x 3 - 4r 2 x 3 - 2r 2 xy 2 
+ 2 r 4 x + 3 r 4 xO - 3 x 2 xO r 2 + 2 x 5 ) 

> g5 : =f actor (reducedG [5] ) ; 

g5 := 8xp y 2 — 6 x 3 p + 7 xOpx 2 + 6 xpr 2 + 4px0 y 2 + y xO x 2 — y xO r 2 + 8yxp 2 
-12yx0p 2 + 2x0pr 2 - 8xp 3 + 8 xO p 3 

> g6 : =f actor (reducedG [6] ) ; 

g6 := p(8p 2 r 2 - 8y 2 p 2 -4x0xp 2 +4x 2 p 2 + 4xyx0p + 8py 3 + 2ypx 2 -8ypr 2 
+ 3 xO x r 2 - 2 r 2 y 2 + 2 x 4 + 2 r 4 + 2 x 2 y 2 - 4 x 2 r 2 - 3 xO x 3 - 4 x y 2 xO) 

> g7 : =f actor (reducedG [7] ) ; 

g7 := xOx 4 -2x 2 xO r 2 +44x 2 p 2 xO + 32 y 2 xO p 2 - YlyxOpr 2 - 80 yxOp 3 + r 4 xO 

+ 16 xOp 2 r 2 + 48 xOp 4 - 2 x 3 yp - 38 x 3 p 2 - 40 x p 2 y 2 + 2 xypr 2 + 56 yxp 3 
+ 20 x p 2 r 2 - 48xp 4 

> g8 : =f actor (reducedG [8] ) ; 

g8 :=3x0 2 r 2 - 12 xO 2 p 2 + 16 x y xO p + 32 xO xp 2 - 20x 2 p 2 - 48 y 2 p 2 + 48p 2 r 2 

- 2 xO x 3 + 2 xO x r 2 + 4 y p x 2 

> g9 : =f actor (reducedG [9] ) ; 

g9 :— -4py 2 - 4px 2 + 6pxx0 - 2px0 2 + 4pr 2 + xO 2 y 

8 We display these polynomials in a form returned by Maple. Polynomials gi, 92, 93i 94, 96 can be further simplified by dividing 
them by p under the assumption that p ^ 0. 
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> glO :=f actor (reducedG [10] ) ; 

glO := -2 xO x 2 + 3 xO 2 x + 2 xO r 2 + 4 y xp - 8px0 y - 8 xp 2 + 8 xO p 2 

> gll:=factor(reducedG[ll] ) ; 

gll := xO 3 - 8xp 2 + 8 xO p 2 - 4px0y 

> gl2 : =f actor (reducedG [12] ) ; 

gl2 := ApyO - xO 2 

> gl3 : =f actor (reducedG [13] ) ; 

glS :=-2xp + 2x0p- xOy + xO yO 

> gl4 : =f actor (reducedG [14] ) ; 

gU := y 2 -2yy0 + yO 2 + x 2 - 2 x xO + xO 2 - r 2 

B Appendix: Reduced Grobner basis for the ideal 

I = (fh /2, /3> C R[y , x , x, y, a, b, r] 

For the ideal 

I = {b 2 yl+a 2 xl-a 2 b 2 ,y 2 -2yy a +yl+x 2 -2xx Q +xl-r 2 , b 2 y x-b 2 y x -a 2 xoy+a 2 xoyo) C R[yo, x , x, y, a, b, r] 

using lex order with y >x >x>y>a>b>r, a reduced Grobner basis for I is given by the following 
fifteen polynomials: 9 

> gl :=f actor (reducedG [1] ) ; 

9 We display these polynomials in a form returned by Maple. All polynomials except gis and gi4 can be further simplified 
by dividing them by a 2 , b 2 under the assumption that a, b ^ 0. 
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gl :=a 2 b 2 (6 b 4 x 4 a 4 - 4 y 6 b 4 a 2 - 6 a 6 x 4 b 2 - 4 6 6 x 2 a 4 + 6 b 4 x 2 a 6 + a 4 b 8 - 2 y 2 o 8 a 2 
+ 6 y 2 & 6 x 2 a 2 + 6 a 6 y 2 x 2 b 2 -10b 4 x 2 a 4 y 2 + y 8 b 4 + b 8 r 4 + b 4 r 8 - 2 6 6 r 6 

- 2 b 4 y 2 x 4 r 2 - 6 x 2 y 4 b 4 r 2 + 6 a; 2 o 4 r 4 y 2 - Ab 4 y 6 r 2 + 2b 4 y 6 x 2 + y 4 b 4 x 4 

+ b 4 x 4 r 4 -2 b 4 r 6 x 2 - 2 b 6 x 2 r 4 + 6 y 4 b 4 r 4 - 2 b 8 y 2 r 2 + 6 6 6 r 4 y 2 - 4 & 4 r 6 y 2 
+ 6 a 4 y 2 o 6 + 6 a 6 y 2 fe 2 r 2 + 2 a 4 fe 2 x 4 y 2 - 4 a 4 x 6 fo 2 + 4 a 6 x 2 6 2 r 2 - 8 6 4 x 2 a 4 r 
+ 10 a 4 x 4 6 2 r 2 — 6 o 4 r 2 a 2 x 2 y 2 — 6 a 4 x 2 b 2 r 2 y 2 — 8 a 4 x 2 6 2 r 4 — 6 o 4 r 2 a 2 x 4 
+ 4fe 4 r 4 a 2 x 2 -6y 4 a 4 b 2 x 2 -6b 4 x 4 a 2 y 2 + 6 6 6 x 2 r 2 a 2 + 2b 4 x 2 y 4 a 2 
+ 10 b 4 y 4 a 2 r 2 + 4 b 6 y 2 r 2 a 2 - 8 b 4 y 2 r 4 a 2 - 8 b 4 y 2 a 4 r 2 + 2 r 6 a 2 fe 4 + a 8 6 4 

- 6 a 4 y 4 r 2 b 2 + 4 a 4 y 2 r 4 b 2 - 2 a 8 b 2 x 2 + x 4 a 4 y 4 + 4 x 2 b 6 y 2 r 2 - 6 y 4 b e a 2 
-4a 6 y 2 b 4 -2y 4 a 4 x 2 r 2 + 2 b 2 r 2 a 2 x 4 y 2 -2r 8 a 2 b 2 + r 8 a 4 + a 8 x 4 + a 4 x 8 

- 4 a 4 x 6 r 2 + 6 a 4 x 4 r 4 + 2 fe 2 r 2 a 2 x 6 - 6 fe 2 r 4 a 2 x 4 - 2 a 6 r 6 + 2 a 6 x 6 + a 8 r 4 
+ 2 y 6 b 2 a 2 x 2 - 2 r 6 a 4 y 2 + 6 r 6 6 2 y 2 a 2 - 2 a 8 x 2 r 2 + y 4 a 4 r 4 - 6 y 4 b 2 r 4 a 2 
-2a 6 y 2 r 4 + Ax 2 a 6 y 2 r 2 + 6 x 2 a 4 y 2 r 4 + 2 y 4 b 2 x 2 a 2 r 2 - 10 x 2 b 2 y 2 r 4 a 2 

+ 2yH 2 a 2 r 2 + 6 r 6 b 2 x 2 a 2 + 4 b 2 x 4 y 4 a 2 + 2 a 4 x 6 y 2 - 6 a 6 x 4 r 2 - 4 r 6 a 4 x 2 
+ 2 6 2 x 6 a 2 y 2 + 6 r 4 a 6 x 2 - 6 a 4 x 4 r 2 y 2 + 2 6 6 y 6 - 2 a 8 6 2 r 2 - 2 x 2 y 4 b 6 
+ 6 a 4 y 4 b 4 + y 4 b 8 - 6 y 4 & 6 r 2 - 2 x 4 a 6 y 2 - 2 a 6 6 6 + 2 6 6 r 4 a 2 - 2 6 8 r 2 a 2 
+ 2 b 4 a 6 r 2 - 6 6 4 a 4 r 4 + 2 r 2 a 4 fe 6 + 2 r 4 a 6 6 2 + 2 r 6 a 4 b 2 ) 

> g2 : =f actor (reducedG [2] ) ; 

g2 := b 2 (-2ya 2 b 2 + b 2 yO y 2 + yO a 2 b 2 - b 2 yO r 2 + 3 a 2 y x xO - a 2 yO x 2 + a 2 y 2 yO 

- a 4 yO + a 2 yO r 2 + a 4 y - a 2 y 3 — a 2 y x 2 + a 2 yr 2 ) 

> g3 : =f actor (reducedG [3] ) ; 

g3 := a 2 (-x0 b 2 y 2 + xO a 2 x 2 + xO a 2 b 2 - xO a 2 r 2 + 3 y fe 2 yO x - 2 x a 2 b 2 + b 2 x 2 xO 
-b 4 x0 + xO b 2 r 2 +b 4 x-b 2 xy 2 - b 2 x 3 + x6 2 r 2 ) 

> g4 : =f actor (reducedG [4] ) ; 

g4 := 6 2 (3 xO a 2 x 3 + 6 a 4 x xO - 3 b 2 x a 2 xO - 3 x xO a 2 r 2 + 4 a 2 y yO x 2 + 2 a 2 y 3 yO 

-2a 4 yy0 - b 2 yO y 3 + 5 y yO a 2 b 2 - 4 a 2 y yO r 2 + y b 2 yO r 2 - 3 x 4 a 2 - 5 x 2 y 2 
~2y 4 a 2 -3a 4 x 2 + 2a 4 y 2 + 36 2 x 2 a 2 - y 2 a 2 b 2 -3a 4 b 2 + 6x 2 a 2 r 2 + 5y 2 a 2 
+ 3a 4 r 2 +3b 2 a 2 r 2 -3r 4 a 2 ) 

> g5 : =f actor (reducedG [5] ) ; 

g5 := b 2 (—b 4 x yO r 2 + a 2 b 4 yO x + 3 xO b 4 y a 2 — 5 b 4 x y a 2 + b 4 y 2 yO x + 2 b 2 x 3 a 2 y 

- 7 b 2 x a 2 y 2 yO - 3 xO b 2 y a 2 r 2 - a 2 b 2 yO x 3 - 3 xO a 4 y b 2 - 2 b 2 x a 2 y r 2 
+ 2 b 2 x a 2 y 3 + 5 x a 4 y b 2 + 3 xO b 2 y 3 a 2 + a 4 x y 2 yO - x 3 y a 4 + a 4 x yO r 2 
+ 3 a 4 xO y r 2 — a 6 x yO — a 4 x y 3 — a 4 yO x 3 + y a 4 x r 2 + a 6 x y) 

> g6 : =f actor (reducedG [6] ) ; 

g6 := b 2 (-2 a 4 y 4 - 6 b 2 r 2 a 2 x xO + b 4 y 2 a 2 + 4 a 4 y yO b 2 + 4 yO b 4 y a 2 - b 2 x 2 a 2 y 2 
+ 3b 2 x 2 a 2 r 2 +4 b 2 y 2 a 2 r 2 -Aa 2 y b 2 yOx 2 -4a 2 y b 2 yO r 2 - 2 b 4 yO y 3 
+ 3 b 4 r 2 a 2 - 3 b 2 r 4 a 2 - 3 a 4 b 4 + 2 y b 4 yO r 2 - b 2 y 4 a 2 + 6 a 4 b 2 x 2 + 6 a 4 b 2 r 2 

- 3 a 6 fe 2 + 4 a 4 y yO x 2 - Aa 4 yy0r 2 + 3 a 6 r 2 - 3 a 4 r 4 + 2 a 4 y 3 yO -2a 6 yy0 
+ 6 a 4 x 2 r 2 + 5 a 4 y 2 r 2 + 6 a 6 x xO - 3 x 4 a 4 - 3 x 2 a 6 - 5 a 4 x 2 y 2 + 2 a 6 y 2 ) 
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> g7:=factor(reducedG[7] ) ; 

gl := b 2 (-21 a 4 x b 4 - 6 x 3 b 4 a 2 + 2 a 6 y 2 x - 3 x 5 a 4 + 6 b e x a 2 + 9 a 6 6 2 a; + 6 a 4 y 2 6 2 a; 

- 5 6 4 y 2 a; a 2 - 5 y 2 a 4 x 3 - 2 y 4 a 4 a; + 12 a 4 fe 2 x 3 - 3 a 6 x 3 - 6 a 6 b 2 xO 

- 6b 6 xO a 2 + 4y a 4 yO x 3 + 2y 3 a 4 x yO - 2y a 6 x yO - 14y a 4 b 2 yO x 

- Ay a 4 xyO r 2 + 5 y 2 a 4 a; r 2 - 4 yOb 2 x 3 a 2 y + 22 yO b 4 xy a 2 + 6 a 4 a; 3 r 2 
+ 3 a 6 x r 2 - 2 6 4 x yO y 3 - 3 a 4 x r 4 - 6 2 a; 3 a 2 y 2 - 6 2 y 4 x a 2 - 6 b 4 y 2 a 2 xO 

- 12 b 2 a 4 xtf r 2 + 6 b 4 x 2 a 2 xO + 6 xtf 6 4 r 2 a 2 + 6 xtf a 6 r 2 + 6 xtf a 4 6 2 y 2 
+ 2b 4 xyy0r 2 + 12x0 a 4 6 4 + 3 b 2 x 3 a 2 r 2 - Ay a 2 b 2 x yO r 2 -3b 2 xr 4 a 2 
-6b 2 x 2 r 2 a 2 x0 + 4 b 2 y 2 xa 2 r 2 + 9 6 4 r 2 a 2 x) 

> g8:=factor(reducedG[8] ) ; 

g8 := fe 2 (4 a 4 y 3 - 2 a 4 x 2 y - 4 a 2 y 3 x 2 - 2 y 5 a 2 - 2 a 6 y - 2 a 2 b 2 yO x 2 + 4 a 2 b 2 yO y 2 

+ 4a 2 b 2 y0r 2 + 2a A yb 2 - 2b 2 y 3 a 2 -2b A ya 2 - b A yO r 2 + b 4 yO a 2 + b 4 yO y 2 

-3a 4 y0b 2 + 2a 4 yr 2 + 2b 2 ya 2 r 2 +4b 2 ya 2 x 2 + 2a 6 y0 + 4x 2 a 2 yr 2 

+ 3 x 2 a 2 y 2 yO - 2 x 2 a 2 yO r 2 - 3 a 2 y 2 yO r 2 - b 2 y 2 yOx 2 + 2 b 2 y 2 yO r 2 

+ b 2 x 2 yO r 2 + 4 y 3 a 2 r 2 - 2 y r 4 a 2 + r 4 a 2 yO - 3 a 4 yO r 2 - b 2 y 4 yO - b 2 r 4 yO 

+ 2 a 2 y 4 yO + a 2 yO x 4 + 3 x 2 a 4 yO - 4 a 4 y 2 yO -2x 4 ya 2 ) 

> g9:=factor(reducedG[9] ) ; 

g9 := b 2 (-4 b 2 x 2 a 2 y 2 yO - 2 x 2 a 2 b 2 yO r 2 - 2 b 4 y 4 yO + b 4 y 3 a 2 - b 2 y 5 a 2 -3ya 4 b 4 
+ 4b 4 y 2 yOa 2 +4 b 4 y 2 yO r 2 - b 2 x 2 a 2 y 3 + 6 x 2 a 4 y b 2 - 2 a 2 y 2 b 2 yO r 2 
+ 2a 4 y 2 yO b 2 + 2 b 2 y 3 a 2 r 2 - yb 4 r 2 a 2 - yb 2 r 4 a 2 +8ya 4 b 2 r 2 + b 2 x 2 a 2 yr 2 
+ 2 r 4 a 2 b 2 yO + 2 a 2 b 4 yOr 2 -2 a 6 yOb 2 + a 6 yb 2 -2 r 4 b 4 yO - 3 a 4 x 4 y 
+ 4a 6 y 3 - 2 a 4 y 5 - 5 x 2 a 4 y 3 - 2 a 6 r 2 yO + a 6 r 2 y + 5 a 4 r 2 y 3 - 3 a 4 r 4 y 
+ 2a 4 y 4 y0 -4a 6 y 2 y0 + 2 a 6 x 2 yO + 2 a 8 yO + 6 a 4 x 2 y r 2 + 4 a 4 x 2 y 2 yO 

- 4 a 4 r 2 yO y 2 - a 6 x 2 y - 2 a 8 y) 

> glO :=f actor (reducedG [10] ) ; 

glO := b 2 (2 a 6 b 4 - 2 a 4 b 6 - 16 a 4 y 4 b 2 + 2 a 4 x 6 - 4 b 4 r 4 a 2 + 2 b 6 r 2 a 2 - 4 b 6 yO y 3 

+ 3 a 4 x 4 y 2 - 2 a 8 y 2 + 4 a 6 y 4 - 4 a 6 r 2 b 2 + 4 a 4 b 2 x 2 y 2 + 9 a 6 y 2 b 2 + 2b 4 a 4 r 2 
+ 6 b 4 x 2 a 4 - 6 y 2 a 4 x 2 r 2 + 2 r 6 6 2 a 2 - 4 x 2 a 6 r 2 + 6 x 2 a 4 r 4 + 6 6 2 x 4 y 2 a 2 
-6 a 4 x 4 r 2 +3y 2 a 2 r 4 b 2 + 13 y 4 b 2 a 2 x 2 - 12 y 4 b 2 a 2 r 2 - 2 r 6 a 4 + 2 r 4 a 6 
+ a 4 y 4 r 2 + 3a 4 y 2 r 4 + 5 a 6 y 2 r 2 + 2a 6 x 4 + 2 b 2 r 2 a 2 x 4 - 4 x 2 b 2 r 4 a 2 
-9y 2 b 2 r 2 a 2 x 2 + 7 y 6 a 2 b 2 - 6 x 4 a 4 b 2 - 4 b 4 x 2 r 2 a 2 + 2 y 5 a 4 yO - 4 y 3 a 6 yO 
+ 2ya 8 y0 + 6b 4 y 5 y0 + 6b 2 r 2 a 2 x 2 yy0-4: b 2 x 2 a 2 y 3 yO + 14 6 2 y 3 a 2 yO r 2 

- 6b 2 a 2 yr 4 yO - 4b 4 x 2 y yO r 2 + 4 6 4 x 2 y 3 yO - 8 b 2 y 5 a 2 y0 + 6 a 6 x 2 y yO 

- 6 a 6 y r 2 yO - 12 b 4 y 3 yO r 2 + 6 b 4 y r 4 yO + 2 a 4 b 2 r 4 - 2 y 6 a 4 + 18 a 4 y 3 b 2 yO 

- 10 a 6 y y0b 2 + 8 a 4 yO b 4 y + b 4 y 2 a 2 r 2 + 4b 6 y yO r 2 - 8 a 4 y b 2 yO x 2 

+ 28 a 4 y b 2 yO r 2 - 26 a 2 y b 4 yO r 2 + 4 a 4 b 2 x 2 r 2 - 18 a 4 b 2 y 2 r 2 - 12 a 2 fe 4 yO y 3 

- x 2 y 4 a 4 - 4 x 2 a 6 6 2 + 3 y 4 6 4 a 2 - 7 a 4 6 4 y 2 + 6 y 2 fe 6 a 2 - 12 b 4 x 2 y 2 a 2 

- 5a 6 y 2 x 2 ) 

> gll:=factor(reducedG[ll] ) ; 
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gll := b 2 {-12 b 4 x 4 a 4 + 18 y 6 b 4 a 2 + 14 6 6 x 2 a 4 - 6 b 4 x 2 a 6 - 5 a 4 b s + 15 y 2 b s a 2 

- 27 y 2 b 6 x 2 a 2 + 6 a 6 y 4 b 2 ~8 a 6 y 2 x 2 b 2 - 6 b 4 x 2 a 4 y 2 - 11 a 4 y 2 b 6 

+ 22 a 6 y 2 b 2 r 2 + 15 a 4 6 2 x 4 y 2 + 24 6 6 y 5 yO + 2 a 4 x 6 b 2 + 16 a 6 x 2 b 2 y yO 

- 24 b 4 x 2 a 4 yy0~8 a 6 x 2 6 2 r 2 + 6 & 4 x 2 a 4 r 2 ~ 9 a 4 x 4 6 2 r 2 - 33 b 4 r 2 a 2 x 2 y 2 
+ 18 b 4 r 2 a 2 x 2 y yO - 29 a 4 x 2 b 2 r 2 y 2 - 4 a 4 x 2 fc 2 r 2 y ytf + 12 a 4 x 2 b 2 r 4 

+ 3b 4 r 2 a 2 x 4 - 10 b 4 r 4 a 2 x 2 - 6 6 4 x 2 a 2 y 3 yO + 8 x 2 b 6 y 3 yO - 8 x 2 6 6 y yO r 2 

- 42 b 6 y 3 yO a 2 - 48 b e y 3 yO r 2 - 7 y 4 a 4 b 2 x 2 + 50 6 4 y 3 a 2 yO r 2 + 9 6 4 x 4 a 2 y 2 

- 9 fe 6 x 2 r 2 a 2 + 43 b 4 x 2 y 4 a 2 - 22 b 4 y 5 a 2 yO + 50 fe 4 y 3 a 4 yO - 29 b 4 y 4 a 2 r 2 

+ 10 b 6 y 2 r 2 a 2 + 4 b 4 y 2 r 4 a 2 - 73 b 4 y 2 a 4 r 2 + 7 r 6 a 2 b 4 + a 8 b 4 + 86 y a 4 fe 4 yO r 2 

- 78ya 2 6 6 yO r 2 - 10y 3 6 8 yO - 28ya 6 6 4 yO + 24ya 4 b 6 yO + 10y& 8 yO r 2 

- 28 a 2 y r 4 b 4 yO + 4 a 4 y 3 r 2 b 2 yO - 18 a 6 y r 2 yO 6 2 + 4 a 4 y r 4 6 2 yO 

- 9 a 4 y 4 r 2 b 2 + 14 a 4 y 2 r 4 b 2 + 24 y fe 6 r 4 yO - 2 a 8 fe 2 x 2 - 6 y 3 a 6 yO 6 2 

+ 4ya s y0 b 2 + x 4 a 4 y 4 + y 4 b 6 a 2 - 4a 8 y 2 6 2 + 22 a 6 y 2 b 4 -2y 3 b 2 r 2 a 2 yO x 2 

- 2 y 4 a 4 x 2 r 2 - 2 x 4 b 4 yyOr 2 - b 2 r 2 a 2 x 4 y 2 - r 8 a 2 b 2 + r 8 a 4 + a 8 x 4 - 2 y 8 6 2 a 2 
+ a 4 x 8 - 4 a 4 x 6 r 2 + 6 a 4 x 4 r 4 + 6 2 r 2 a 2 x 6 - 3 6 2 r 4 a 2 x 4 - 2 a 6 r 6 + 2 a 6 x 6 + a 8 r 4 

- y 6 6 2 a 2 x 2 + 2 x 2 y 6 4 r 4 yO - 2 r 6 a 4 y 2 + 5 r 6 6 2 y 2 a 2 - 2 a 8 x 2 r 2 + y 4 a 4 r 4 
-9y 4 b 2 r 4 a 2 + 2 a 2 y 7 b 2 y0-2a 6 y 2 r 4 - 4y 5 b 2 a 2 yO r 2 + 2 y 3 a 2 r 4 b 2 yO 

+ 4x 2 a 6 y 2 r 2 +6 x 2 a 4 y 2 r 4 + 5 y 4 6 2 x 2 a 2 r 2 - 7 x 2 6 2 y 2 r 4 a 2 + 2 y 5 fe 2 a 2 yO x 2 

- 4 x 2 b 4 y 3 y0r 2 + 7 y 6 b 2 a 2 r 2 + 3 r 6 6 2 x 2 a 2 + 2 x 2 6 4 y 5 yO + 4 6 2 x 4 y 4 a 2 
+ 2a 4 x 6 y 2 -6 a 6 x 4 r 2 -4r 6 a 4 x 2 + 2x 4 fo 4 y 3 y0 + 3b 2 x 6 a 2 y 2 + 6r 4 a 6 x 2 

- 6 a 4 x 4 r 2 y 2 - 2 a 8 b 2 r 2 - 41 a 4 y 4 b 4 - 2 x 4 a 6 y 2 + 4 a 6 fe 6 - 11 fe 6 r 4 a 2 
+ 5 b s r 2 a 2 - 10 b 4 a 6 r 2 + 2 fe 4 a 4 r 4 + 7 r 2 a 4 b 6 + 8 r 4 aH 2 - 5 r 6 a 4 fe 2 ) 

> gl2 :=f actor (reducedG [12] ) ; 
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b 2 (-3 y 6 b 4 a 2 x 2 - 9 a 4 r 6 o 4 + 2 a 6 6 8 + 2 a 8 r 4 o 2 - 19 a 6 r 2 b 6 + 32 a 6 fe 4 r 4 

- 13 a 8 o 4 r 2 + 20 a 4 6 s r 2 - 10 a 4 o 6 r 4 + 8 a 8 o 6 + 24 y 2 b w a 2 - 15 a 6 r 6 b 2 
+ 10 a 6 x 6 b 2 - 17 a 6 x 4 y 2 b 2 - 39y 2 b s x 2 a 2 + 23 a 6 y 2 x 2 b 4 - 16 a 8 x 2 y 2 b 2 

- 37 a 4 y 4 x 2 o 4 + 24 x 2 y 4 a 6 6 2 - 16 b w y 3 y0 + 2 b 4 x 6 a 4 + 94 y 4 6 6 x 2 a 2 

+ 9y 2 b 6 x 4 a 2 - 52 o 6 x 2 a 4 y 2 + 48 o 4 x 4 a 4 y 2 - 8 b 10 a 4 - 18 a 6 o 4 x 4 + 8 b 10 r 2 a 2 

- 22 b s r 4 a 2 + 21 b 6 r 6 a 2 + 16 b w y yO r 2 - 8 r 8 a 2 o 4 + 8 r 8 a 4 6 2 - 16 b 6 x 4 a 4 
+ 21 b 8 x 2 a 4 + 2 y 3 o 4 r 2 a 2 yO x 2 - 20 y 4 a 4 o 2 x 2 r 2 - 2 a 10 6 4 - 2 a 10 x 4 

+ 4 a 10 b 2 x 2 + a 4 x 10 - 3 a 8 x 6 - 19 a 8 6 4 x 2 + 14 a 8 x 4 fe 2 - 16 a 4 o 2 y 6 x 2 

- 2 a 6 x 6 y 2 + 5 y 4 a 4 x 6 + 4 a 4 x 8 y 2 - 6 x 4 a 6 y 4 + 6 x 4 a 8 y 2 + 2 y 6 a 4 x 4 + 6 6 6 x 2 

- 14 y 4 o 8 a 2 + 24 6 6 y 6 a 2 + 20a 6 y 2 b 6 -4x 4 o 6 y ytf r 2 - 8 6 4 r 2 a 2 x 4 y 2 

- 8 a 4 x 6 o 2 r 2 + 24a 4 x 4 b 2 r 4 + 6 4 r 2 a 2 x 6 - 10 6 4 r 4 a 2 x 4 + 4x 4 6 6 y 3 ytf 

+ 3 6 4 x 6 a 2 y 2 + 4 b 4 a 4 x 2 r 4 + 40 r 4 a 6 x 2 6 2 - 10 r 2 a 4 x 2 b 6 + 17 r 6 o 4 x 2 a 2 

- 24 r 6 a 4 x 2 o 2 - 20 a 4 x 4 6 2 r 2 y 2 - 30 x 2 b 4 y 2 r 4 a 2 - 2 x 2 a 2 y r 4 b 4 yO 

+ 40 x 2 a 4 y 2 r 4 o 2 - 12 b 6 x 2 r 4 a 2 + 8 x 2 b 6 y 5 yO + 6 b 4 x 4 y 4 a 2 + 3 6 6 x 4 a 2 r 2 

- 35 a 6 x 4 b 2 r 2 + 16 y 4 o 4 x 2 a 2 r 2 + 1 16 y 3 a 2 6 6 yO r 2 - 4 y 5 6 4 a 2 ytf r 2 
+ 2 y 3 a 2 r 4 o 4 yO + 34 x 2 a 6 y 2 o 2 r 2 - 20 x 2 b 6 y 3 yO r 2 + 10 6 8 x 2 yO y 3 

- 13o 8 x 2 r 2 a 2 - 50 y 5 b 6 yO a 2 - 10 y 5 b 6 yO r 2 +39y 6 6 4 a 2 r 2 - 8y 4 a 4 b 2 x 4 

- 67 y 4 a 4 o 4 r 2 - 31 y 4 o 6 a 2 r 2 - 8 y 3 a 6 b 4 yO + 8 y 3 o 6 r 4 yO + 24 y 4 a 6 r 2 o 2 

- 16 y 6 a 4 r 2 o 2 + 28 y 4 a 4 r 4 o 2 - 56 y 4 o 4 r 4 a 2 + 3 b 4 a 4 x 4 r 2 + 6 o 4 a 6 x 2 r 2 
+ 4 a 8 y 2 6 4 - 10 y 8 6 4 a 2 + 100 a 4 y 3 o 6 ytf - 17 a 6 y 2 6 2 r 4 - 16 a 8 y 2 6 2 r 2 

- 178 a 4 o 6 y 2 r 2 + 48 a 4 6 4 r 4 y 2 - 94 6 s y 3 yO a 2 - 114 o 8 y 3 ytf r 2 + 32 b s y 2 a 2 r 2 

- 14 b 6 r 4 y 2 a 2 - 20 r 6 a 4 o 2 y 2 + 35r 6 b 4 y 2 a 2 + 6ya s b 4 y0 - My a 6 b 6 yO 
+ 56 y 6 8 r 4 yO - 12 6 4 y a 4 x 2 yO r 2 + 30 6 4 y a 6 x 2 yO - 2 6 6 y 3 a 2 yO x 2 

- 48 6 6 y a 4 yO x 2 + 32 6 6 y a 2 x 2 yO r 2 - 82 6 4 y 2 a 4 x 2 r 2 - 10 6 8 x 2 y y# r 2 

- 70 6 6 y 2 x 2 a 2 r 2 + 8 a 4 y 3 o 4 yO r 2 - 34 y a 6 6 4 yO r 2 + 170 y a 4 o 6 yO r 2 
+ 10 y a 4 r 4 b 4 yO - 66 y a 2 6 6 r 4 ytf - 152 y o 8 a 2 yO r 2 + 77 a 6 y 2 o 4 r 2 

+ 58 6 8 y 5 yO + 48 y a 4 o 8 yO - 52 a 4 y 4 b 6 + 2 6 4 y 9 yO + 4 b 6 y 7 yO + 6 6 4 y 7 ytf x 2 

- 12 a 6 y 4 o 4 + 12 x 2 a 6 y 4 r 2 + 10 x 2 y 3 6 4 r 4 yO - 12 x 2 a 8 y 2 r 2 + 4 r 8 a 4 y 2 

- 7 r 8 6 2 y 2 a 2 + 2 y 6 a 4 r 4 - 2 y 3 6 4 r e yO - 4 y 6 a 4 x 2 r 2 - 6 o 4 y 7 yO r 2 

+ 1 1 b 2 y 6 x 4 a 2 + 4 6 2 y 8 x 2 a 2 + 4 b 2 y 8 a 2 r 2 - IS b 2 y 6 r 4 a 2 + 6 b 4 y 5 r 4 ytf 

+ 5 r 8 a 4 x 2 - 4 r 8 o 2 x 2 a 2 - 16 a 8 6 2 x 2 r 2 + 9 a 8 x 4 r 2 - 9 a 8 r 4 x 2 - 5 a 4 x 8 r 2 

+ 10 a 4 x 6 r 4 + b 2 r 2 a 2 x 8 - 4 6 2 r 4 a 2 x 6 + 2 x 6 6 4 y 3 yO + 3 b 2 x 8 a 2 y 2 

+ 6 r 6 b 2 x 4 a 2 - 10 r 6 a 4 x 4 + 6 x 4 6 4 y 5 yO + 10 6 2 x 6 y 4 a 2 - 12 x 4 b 2 y 2 r 4 a 2 

+ 24 x 4 a 4 y 2 r 4 - 5 y 4 b 2 x 4 a 2 r 2 - 10 x 4 b 4 y 3 yO r 2 - 6 a 6 x 2 y 2 r 4 - 15 y 4 a 4 x 4 r 2 

- 2x 6 b 4 yy0r 2 - 2b 2 r 2 a 2 x 6 y 2 + 4r 4 x 4 b 4 yyO + 15 x 2 y 4 a 4 r 4 

- 20 x 2 y 4 fe 2 r 4 a 2 - 16 x 2 r 6 a 4 y 2 + 18 x 2 r 6 o 2 y 2 a 2 + 12 x 2 y b 6 r 4 yO 

+ 6 a 6 x 4 r 2 y 2 - 16 a 4 x 6 r 2 y 2 - 14 b 4 y 5 x 2 yOr 2 + 2 b 2 y 6 x 2 a 2 r 2 - 6 y 4 a 6 r 4 

- 5 y 4 r 6 a 4 + 15 y 4 r 6 b 2 a 2 + 6 a 8 y 2 r 4 + 2 a 6 r 6 y 2 + 2 y b 4 r 6 a 2 y0-2y 6 6 r 6 yO 
+ 20 y 6 a 4 b 4 + 4 a 10 x 2 r 2 - 2 r 6 x 2 b 4 y yO - 2 a 10 r 4 + 4 a 10 6 2 r 2 + 3 a 8 r 6 

+ r 10 a 2 b 2 - r 10 a 4 ) 
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> gl3 :=f actor (reducedG [13] ) ; 

gl3 :=y 2 - 2yy0 + y0 2 +x 2 -2xx0 + xO 2 - r 2 

> gl4 : =f actor (reducedG [14] ) ; 

gU := b 2 y0x- b 2 yO xO - a 2 xO y + a 2 xO yO 

> gl5 : =f actor (reducedG [15] ) ; 

gl5 := b 2 (xO yO x + y yO 2 - yO x 2 - yO y 2 - yO a 2 + yO r 2 + y a 2 ) 

C Appendix: Reduced Grobner basis for the ideal 

Is = (h h h 2 , h 3l /i 4 , h 5 , g) C R[x, y, a, 6, r] 

In order to solve system (|38p. we have found a reduced Grobner basis for using the total term order tdeg 
for y>x>a>b>r with ties broken by inverse lexicographical order. We used package Groebner from 
Maple 8 supplemented by a small custom package RJGrobner written by the Authors. [2] The reduced basis 
contains fourteen homogeneous polynomials with degrees ranging from 10 through 15 and the number of 
monomials terms ranging from 51 through 112. We will not display these polynomials here. However, here 
is the factorization the fourth polynomial 64 in the basis: 

b 4 = b 2 x(4 y 2 r 4 a 2 - 5 y 4 a 2 r 2 + a 4 r 2 y 2 - 2 b 2 x 2 r 4 + b 2 x 4 r 2 ~ 9 x 2 b 2 y 2 a 2 - 11 b 2 y 2 a 2 r 2 
-2x 2 b 4 r 2 +4a 6 r 2 -2a s -b e a 2 -a 4 b 4 + Aa 6 b 2 + 6 y 2 a 6 - a 4 r 4 - r e a 2 

- Ay 4 b 4 - 6 y 4 a 4 - 4a 6 x 2 + 2 if a 2 - a 4 x 4 - b 6 y 2 + b 6 r 2 + r 6 b 2 - Ab 4 r 4 

+ 2a 4 b 2 x 2 -9a 4 b 2 y 2 + 2 x 2 b 4 y 2 + 5 b 4 y 2 a 2 - a 4 x 2 y 2 - 3 b 2 a 2 x 4 + 6 b 2 y 4 a 2 (79) 

io2 7.4 2 6 1.2 1 26 r»i,2 222ir»422 n242,o224 

+ 3 x a — y b + a x — 2 x a r + 2 a a; r — 3 a x r + 3 a x r 

o 2 2 2 2 1 ,1 1.2 2 2 2,ol.4 2 2 , c l2 4 2,t-i,4 2 2 ol2 4 2 

— 8 x y a r + 4 o x y r +8 b y r + 5 o r a + 5 o r a — 3 o r y 

iOJ.2 42 -inl.2 42 7.2 42 r ) 4j.2 2 l/ |242,r242\ 

+ 3b y r — 10 bar — b x y —2ybx + 4 a x y + 5 a y x ) 

It appears, that since b > 0, one possible solution is x = 0. Upon substituting x — into the list i 7, shown 
in one gets 

[0, -2/ A(2 6 2 y 4 + 6 4 y 2 - 4 6 2 r 2 y 2 - 4 6 2 y 2 a 2 + y 2 a 2 r 2 - b 4 a 2 - b 4 r 2 + 2 b 2 r 4 + 2 a 4 b 2 -a 4 r 2 -r 4 a 2 ), 
- A(-2 b 2 y 4 + r 2 y 4 - b 4 y 2 + 2 b 2 r 2 y 2 + 4 b 2 y 2 a 2 - 3 y 2 a 2 r 2 - 2 r 4 y 2 + b 4 a 2 - b 4 r 2 - 2 a 4 b 2 

+ 2b 2 a 2 r 2 +r 6 + 2 a 4 r 2 -3r 4 a 2 ), 
(y + r - a) (y + r + a) (a - y + r) (-a -y + r) (-y 2 + r 2 + a 2 - 2 b 2 ) A, 

A(-y 4 a 2 + 2 b 2 y 4 + b 4 y 2 - 4 b 2 r 2 y 2 - 2 b 2 y 2 a 2 + 3 y 2 a 2 r 2 + 2 a 4 y 2 + b 4 a 2 - b 4 r 2 - 2 b 2 a 2 r 2 

+ 2& 2 r 4 + 3a 4 r 2 -2r 4 a 2 -a 6 ), 
(y + r - a) (y + r + a) (a - y + r) (—a -y + r) A 2 ] 

(80) 

where 

A = -y 2 b 2 -b 4 ~r 2 a 2 + b 2 a 2 + b 2 r 2 . (81) 
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It can be easily verified that setting y = —r — a or y = r + a in l|8U|l does not yield any critical values for r, 
as the system becomes, respectively, 

F x = [0, -2(r + a) {b 2 + ra) 4 ra, 2 {b 2 + r a) 4 r (r + a), 0, -2 (b 2 + r a) 4 a (r + a), 0], (82) 
F± = [0, 2(r + a)(b 2 +ra) 4 ra, 2 {b 2 + r a) 4 r (r + a), 0, -2 (b 2 + r a) 4 a (r + a), 0] (83) 

asa>6>r>0 and is insolvable. However, setting y = —r + a or y = r — a does yield r cr u = — and makes 
the system solvable in each case: 

F 2 = [0, 2{r-a){-b 2 + ra) 4 ra, 2 {-b 2 + r a) 4 r (r - a), 0, 2 (-b 2 + r a) 4 a {r - a), 0], (84) 

F 3 = [0, -2{r-a){-b 2 + ra) 4 ra, 2(-6 2 + r«) 4 r(r-a), 0, 2 (-b 2 + r a) 4 a {r - a), 0] (85) 

Another possibility for solving the system is to set 

_ 6 4 _ r 2 fl 2 + 6 2 fl 2 + & 2 r 2 
V = 1~2 ( 86 ) 

Thus, 



/ /_ ^4 _ r 2 fl 2 _|_ ^2 fl 2 _|_ ^2 r 2 \ 

h = °» v = ± v p J ( 87 ) 

provide two solutions to the system without any condition on r except that r > 0. These are the two virtual 
singular points found earlier that are always present for a > b > r > 0. Observe also that the critical value 
of r satisfies 

b 2 fb\ 
Tcrit = — = — ] b < b since a > b. 



a 

Another possibility to solve the system is to set 

y 2 = r 2 + a 2 - 2 b 2 

When we substitute the above expressions for y 2 into F, we get 

[0, -6 y (-b 2 + r a) 2 (b 2 +ra) 2 b 2 ,6 {-b 2 + r a) 2 {b 2 + r a) 2 (-6 + r)(b + r), 0, 

6 (-b 2 + r a) 2 (b 2 + r a) 2 {b -a){b + a), (89) 
(a + r — y) (a + r + y) (i — a + y) (—a + r — y) (—b 2 + r a) 2 (b 2 + r a) 2 ] 

The above shows that in this case in order to satisfy the equations of the variety of singular points V(/g) 

,2 

we must also have r = r cr u = — . Hence we get 

fa 2 -b 2 \ c 2 \ b 2 

x = 0,y = ±{ =— when r=— , (90) 

\ a J a J a 

which is the special case of the above singular points discussed above. 

Finally, we have special cases that we are not interested in but we mention them for completeness. They 
result when we set y = in the list F (in addition to x — 0. The first case is when a = b, that is, when the 
ellipse is a circle. The second case is when r = b or r = a. 

Thus, to summarize, we have these two cases of singular points when x = : 

x = 0, y = fr-f r — Sl - 2 - j and r is arbitrary as long as b > r > 0. This yields singular points 

that are always present. 
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2. (x = 0, y = ±(— r + a) = ± (tt)J when r = r cr n — ^- and c 2 = a 2 — & 2 . This yields singular points 
in the special case r — r cr n. 

In order to complete the analysis of all solutions of the system F, we need to consider the case when 
x, y 7^ 0. The fourteen-polynomial reduced Grobner basis found for Is when solving system (|38fl . could be 
simplified by diving each polynomial by x,y,a 2 , or b 2 , as appropriate. Let / C R[a;, y, a, b, r] denote the 
ideal generated by these fourteen simplified polynomials. In an effort to eliminate variable y from some 
polynomials, we computed 10 the Grobner basis for the first elimination ideal I x = I fl R[x,a,b,r] for the 
lexicographic order y>r>a>b>x. The first polynomial in this new thirteen-polynomial homogeneous 
basis for I x is of degree 12. It is the square of the following polynomial g x C R[x, a, b, r] of degree 6 in x : 

g x = (b 6 -3b 4 a 2 + 3a 4 b 2 - a 6 ) x 6 + (6b 6 a 2 + 3a 6 r 2 -6b 2 a 4 r 2 + 36 4 r 2 a 2 - 3 b s - 3 a 4 b 4 ) x 4 

+ (-21 b 4 a 4 r 2 + 3 6 10 - 3 6 8 a 2 + 21 6 6 r 2 a 2 + 3 a 4 6 2 r 4 - 3 r 4 a 6 ) x 2 (91) 
-36 4 a 4 r 4 -6 12 + 36 8 r 2 a 2 + a 6 r 6 

Upon substituting 2 = x 2 this becomes an easily solvable cubic equation for z that has the following solution: 

_ 3 6 2 a 2 r 2 - r 2 a 2 (abr)^ - 3 6 3 (abr)^ ar + b A (a for) 3 _ 5 

(a6r)i (-a 2 + b 2 ) (abr)i (-a 2 + b 2 ) 

where 8 is given in 1451) . Thus, we have derived the x-coordinates of the singular points 1S5, Sq, St and S& 
listed in glj, g2|>, US, and jSJ). 

To find the y coordinates of these four singular points, we proceed in a similar fashion except that we 
eliminate the x variable. Thus, compute the first elimination ideal I y = I n a, b, r] for the lexicographic 
order x>r>a>b>y. The first polynomial in this new six-polynomial homogeneous basis for I y is of 
degree 12. It is the square of the following polynomial g y C a, b, r] of degree 6 in y : 

g v = (-b 6 + 3 6 4 a 2 - 3 a 4 6 2 + a 6 ) y 6 + (3 6 2 a 4 r 2 + 3 b e r 2 - 3 a 4 6 4 + 6 a 6 6 2 - 3 a 8 - 6 6 4 r 2 a 2 ) y 4 

+ (3b i a 2 r 4 -3b 6 r 4 -21b i a 4 r 2 + 3a w + 21a 6 b 2 r 2 -3a 8 b 2 )y 2 + b e r 6 (93) 
-36 4 a 4 r 4 -a 12 + 3a 8 fe 2 r 2 

Upon substituting z — y 2 this becomes an easily solvable cubic equation for z that has the following solution: 

-a 4 (abr)i + 3 (abr)i bra 3 - 3 b 2 a 2 r 2 + b 2 r 2 (abr)i _ (-a 2 + b 2 )e 
Z= (abr)i(a 2 -b 2 ) 2 = (a br) 2 / 3 (-a 2 + b 2 ) 2 ( } 

where e is given in 1461) . Thus, we have derived the y-coordinates of the singular points S^,Sq,S7 and 5s 
listed in gTJ, {HJ, ESJ), and (gg) . 
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10 We used gbasis procedure from the Groebner package with the total order lexdeg([y] , [r,a,b,x]) indicating that the 
variable y needed to be eliminated. 
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